Higher-Loop Corrections to the Infrared Evolution of a Gauge Theory with Fermions 
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We consider a vectorial, asymptotically free gauge theory and analyze the effect of higher-loop 
corrections to the beta function on the evolution of the theory from the ultraviolet to the infrared. 
We study the case in which the theory contains Nf copies of a fermion transforming according to 
the fundamental representation and several higher-dimensional representations of the gauge group. 
We also calculate higher-loop values of the anomalous dimension of the mass, j m of {pip at the 
infrared zero of the beta function. We find that for a given theory, the values of 7™, calculated to 
three- and four-loop order, and evaluated at the infrared zero computed to the same order, tend to 
be somewhat smaller than the value calculated to two-loop order. The results are compared with 
recent lattice simulations. 



I. INTRODUCTION 

In this paper we investigate how higher-loop correc- 
tions to the beta function affect the evolution of a vec- 
torial, asymptotically free gauge theory (in (3 + 1) di- 
mensions, at zero temperature) from the ultraviolet to 
the infrared. We assume that the theory contains a cer- 
tain number, Nf, of massless Dirac fermions ip transform- 
ing according to a representation R of the gauge group. 
We consider cases where R is the fundamental, adjoint, 
and rank-2 symmetric or antisymmetric tensor represen- 
tation. We also study the effect of higher-loop correc- 
tions to the anomalous dimension 7 m of the fermion mass. 
This work yields more complete information on the na- 
ture of the evolution of the theory from the ultraviolet to 
the infrared, in particular, on the determination of the 
infrared zero of the beta function and the scaling behav- 
ior of the tptp operator in the vicinity of this zero. We 
will give a number of results for a general gauge group G 
but will focus on the case G = SU(iV). 

We denote the running gauge coupling of the the- 
ory as g(n), with a(jj,) — g{n) 2 / (47r), where jj, is the 
Euclidean energy /momentum scale (which will often be 
suppressed in the notation). The property that the 
SU(JV) gauge interaction is asymptotically free means 
that lim^oo a(p) = 0, and, since the beta function 
is negative for small a, it follows that, as the en- 
ergy/momentum scale ll decreases from large values, a 
increases. As ll decreases and the theory evolves into 
the infrared, two different types of behavior may occur, 
depending on the fermion content. In a theory with a suf- 
ficiently small number of fermions in small enough rep- 
resentations R, as ll decreases through a scale A, the 
coupling a exceeds a critical value a_R jC r, depending on 
R, for the formation of bilinear fermion condensates, and 
these condensates are produced. This may or may not 
be associated with an infrared (IR) zero of the two-loop 
beta function at a value a = am; if the two-loop (2£) 
beta function does have an IR zero, am,2t, then this 
type of behavior requires that am,2t > aR yCr |H, 0- As 
ll decreases toward A and a increases toward the IR zero 
of the beta function, the increase of a as a function of 



decreasing ll is reduced. This gives rise to an a that is 
of order unity, but varies slowly as a function of ll. This 
behavior is interestingly different from the behavior of 
the gauge coupling in quantum chromodynamics (QCD). 
As the condensates form, the fermions gain dynamical 
masses of order A, so that in the low-energy effective the- 
ory applicable at scales 11 < A, they are integrated out, 
and the further evolution of the theory into the infrared 
is controlled by the Nf — beta function. 

Alternatively, if the theory has a sufficiently large num- 
ber Nf of fermions and/or if these fermions are in a large 
enough representation R (as bounded above by the re- 
quirement of asymptotic freedom), then the IR zero of 
the beta function occurs at a value smaller than au tCr , 
so that as the scale ll decreases from large values, the 
theory evolves into the infrared without ever sponta- 
neously breaking chirally symmetry. In this latter case, 
the IR zero of the beta function is an exact IR fixed point 
(IRFP), approached from below as ll — > 0. More com- 
plicated behavior occurs in theories containing fermions 
in several different representations Q ; here we restrict to 
the case of fermions in a single representation. For a given 
asymptotically free theory that features an IR fixed point 
at a = am, the value of this IRFP decreases as a func- 
tion of Nf. There is thus a critical value of Nf, denoted 
Nf ;Cr , depending on R, at which am decreases below 
an, cr . This value serves as the boundary, as a function 
of Nf, between the interval of nonzero 1 < Nf < Nf >cr 
where the theory evolves into the infrared in a manner 
that involves fermion condensate formation and associ- 
ated spontaneous chiral symmetry breaking (SxSB), and 
the interval Nf :Cr < Nf < Nf >max , where the theory 
evolves into the infrared without this condensate forma- 
tion and chiral symmetry breaking, with Nf tTnax denoting 
the maximal value of Nf consistent with the requirement 
of asymptotic freedom. 

The anomalous dimension 7 m contains important in- 
formation about the scaling behavior of the operator tpip 
for which m is the coefficient, as probed on different mo- 
mentum scales. In a theory with an am ~ 0(1), it 
follows that j m may also be 0(1), which can produce sig- 
nificant enhancement of dynamically generated fermion 
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masses due to the renormalization-group factor 



rj = exp 



7m(a(/x)) 



(1.1) 



In a phase where no dynamical mass is generated, 7 m 
simply describes the scaling behavior of the xpip operator. 

There are several motivations for the study of higher- 
order corrections to this evolution of an asymptotically 
free theory into the infrared. First, the critical coupling, 
Q!_R,cr, is generically of order unity, and hence there is a 
need to have a quantitative assessment of the importance 
of higher-loop corrections to the evolution of the the- 
ory. Second, besides being of fundamental held-theoretic 
interest, a knowledge of this evolution plays an impor- 
tant role in modern technicolor (TC) models with dy- 
namical electroweak symmetry breaking, in which the 
slow running of the coupling associated with an approx- 
imate infrared zero of the beta function provide neces- 
sary enhancement of quark and lepton masses 0, (re- 
cent reviews include [4|-[6j]), and can reduce technicolor 
corrections to precision electroweak quantities @, H| . In 
addition to the fundamental representation, fermions in 
higher-dimensional representations have been studied in 
the context of technicolor [f| H|. Fermions in higher- 
dimensional representations of chiral gauge groups have 
long played a valuable role in studies of extended tech- 
nicolor (ETC) models that were reasonably ultraviolet- 
complete and explicitly worked out the details of the se- 
quential breaking of the ETC chiral gauge symmetries 
down to the TC group 10]. Recently, there has been 
a considerable amount of effort devoted to lattice stud- 
ies of gauge coupling evolution and condensate forma- 
tion in vectorial SU(iV) gauge theories as a function of 
Nj^jov fermions in both the fundamental representation 
(8|,[l]J-[lj|and higher representations [ll|-[23[ (a recent 
review is [24]). Thus, another important motivation for 
the present work is to provide higher-order calculations 
that can be compared with these lattice studies. 



II. GENERAL THEORETICAL FRAMEWORK 

A. Beta Function 

The beta function of the theory is denoted f3 = dg/dt, 
where dt = dln/i. In terms of the variable 
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(2.1) 



the beta function can be written equivalently as /3 a = 
da/dt, expressed as a series 

.00 oo 

^ = -2aJ2 bt a e = -2a £ b e a e , (2.2) 

l=x e=i 

where I denotes the number of loops involved in the cal- 
culation of bi and b( — bi/(A-K) 1 . Although this series and 



series for other quantities in quantum field theories do not 
have finite radii of convergence but are only asymptotic, 
experience shows that in situations where the effective ex- 
pansion parameter (here, (ot/n) times appropriate group 
invariants) is not too large, the first few terms can pro- 
vide both qualitative and quantitative insight into the 
physics. The first two coefficients in the expansion 
which are scheme- independent, are [25j 



b 1 = ±(llC A -4T f N J ) 



and \2t 



62 = - [MC\ - A{bC A + 3Cf)T f N f ] 



(2.3) 



(2.4) 



Here Cf = C%(R) is the quadratic Casimir invariant for 
the representation R to which the Nf fermions belong, 
Ca = 02(G) is the quadratic Casimiar invariant for the 
adjoint representation, and T/ = T(R) is the trace in- 
variant for the fermion representation R. Higher-order 
coefficients, which are scheme-dependent (27} , have been 
calculated up to four- loop order (13, 29]. Some further 
details are given in Appendix I. 



B. Anomalous Dimension of the tpip Operator 

The anomalous dimension 7 m for the fermion bilinear 
ipip describes the scaling properties of this operator and 
can be expressed as a series in a or equivalently, a: 
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(2.5) 



where c.£ = ct/(4w) is the ^-loop series coefficient. Via 
Eq. (jl.ip . the anomalous dimension j m governs the run- 
ning of a dynamically generated fermion mass. The co- 
efficients ci have been calculated to four-loop order [30I ]. 
The first two are 



ci = 6C/ 



(2.6) 



and 



c 2 = 2C 



(2.7) 



12 1 6 3 J 

For reference, the coefficient C3 is listed in Appendix I. 
Since as Nf approaches Nf max from below, b\ — \ with 
nonzero 62 and hence our — ¥ 0, and since the perturba- 
tive calculations expresses "f m in a power series in a, it 
follows that as ~/ m — > as Nf approaches Nf^ max from 
below. We note that a conjectured beta function that 
directly relates (3 to 7 has been proposed (3lj |. 



III. PROPERTIES OF BETA FUNCTION 
COEFFICIENTS AND APPLICATION TO 
FUNDAMENTAL REPRESENTATION 

In this section we discuss some general properties of 
the beta function coefficients as functions of Nf, and 
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give particular results for the case of fermions in the fun- 
damental representation. In later sections we consider 
fermions in two-index representations. 
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Since we restrict our considerations to an asymptoti- 
cally free theory, we require that, with our sign conven- 
tions, b\ > 0. This, in turn, implies that 



where 



Nf < Nf t max , 

nc A 



Nf. 



AT 
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(3.1) 



(3.2) 



Thus, for fermions in the fundamental representation, 

Nf^rnaxjund = (U/2)N. 



B. 62 and Condition for Infrared Zero of /3 

We next proceed to characterize the behavior of the 
higher-loop coefficients of the beta function, bi with 
I = 2,3,4, and the resultant zero(s) of the beta func- 
tion, in terms of their dependence on Nf. The two- loop 
results are well-known and are included here so that the 
discussion will be self-contained. Since only the first two 
coefficients of the beta function are scheme-independent, 
it follows that, to the extent that one is in a momentum 
regime where one can reliably use the perturbative beta 
function, the zeros obtained from these first two coeffi- 
cients should be sufficient to characterize the physics at 
least qualitatively. When one includes higher-loop con- 
tributions to the beta function, one expects shifts of ze- 
ros, and there are, indeed, generically substantial shifts 
if zeros of the two-loop beta function occur at a ~ 0(1). 
However, if inclusion of three- and/or higher-loop contri- 
butions to j3 leads to a qualitative change in behavior, 
relative to the behavior obtained from the two-loop (3 
function, then the results cannot be considered fully reli- 
able, since they are scheme-dependent. For example, for 
a given gauge group G and fermion content, if the two- 
loop beta function does not have an infrared zero but 
the three-loop beta function does, one could not con- 
clude reliably that this is a physical prediction of the 
theory. Moreover, it should be noted that even if there is 
no zero of the two-loop beta function away from the ori- 
gin, i.e., a perturbative IRFP, the beta function may ex- 
hibit a nonperturbative slowing of the running associated 
with the fact that at energy scales below the confinement 
scale, the physics is not accurately described in terms of 
the Lajjrangian degrees of freedom (fermions and gluons) 

Another general comment is that the expression of the 
beta function in Eq. (I2.2[) is semiperturbative and does 
not incorporate certain nonperturbative properties of the 



physics, such as instantons, whose contributions involve 
essential zeros of the form exp(— kit /a), where K is a nu- 
merical constant. These instanton effects are absent to 
any order of the perturbative expansion in Eq. (|2.2p but 
play an important role in the theory. For example, they 
break the global U(l)^ symmetry [35| and also enhance 
spontaneous chiral symmetry breaking [37|-[39j. Esti- 
mates of the effects of instantons on the running of a in 
quantum chromodynamics have found that they increase 
this running, i.e., they make /? more negative in the re- 
gion of small to moderate a values (36|. If one were to 
model the effect of instantons crudely via a modification 
of j3 such as 



_da 
Pa ~ dt 



-2a' 
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^6,a £ - 1 +Aexp(-^) 



(3.3) 



then, since A > 0, this would have the effect of increasing 
the value of the smallest (nonzero, positive) IR zero am 
of j3. For a given minimal value of a cr .n for condensate 
formation and spontaneous chiral symmetry breaking, 
since at least at the perturbative level am is a decreas- 
ing function of Nf, it would follow that incorporating 
instanton effects would increase the value of Nf cr , i.e., 
would increase the interval in Nf where there is SxSB. 
Furthermore, since instantons enhance chiral symmetry 
breaking, they would tend to reduce the value of a cr ^R, 
which also has the same effect of increasing Nf <cr . We 
shall comment below on how, although the semipertur- 
bative one-gluon exchange approximation to the Dyson- 
Schwinger (DS) equation does not directly include effects 
of confinement or instantons, it may nevertheless yield an 
approximately correct value of Nf tCr because of another 
approximation involved that has the opposite effect on 
the estimate. 

If one knows the beta function calculated to a maximal 
loop order i max , then the equation for the zeros of the 
beta function, aside from the zero at a = 0, is 



1=1 



1=1 



= . (3.4) 



As is clear from Eq. (|3.4|) . the zeros of j3 away from 
the origin depend only on the £ max — 1 ratios 6^/61 for 

2 < P < f 

The coefficients b\ and 62 are linear functions of Nf, 
while 63 and 64 are, respectively, quadratic and cubic 
functions of Nf. With our sign convention in which an 
overall minus sign is extracted in Eq. (|2.2p , each of these 
coefficients is positive for Nf = 0. The coefficients b\ and 
62 are both monotonically and linearly decreasing func- 
tion of Nf . As Nf increases sufficiently, 62 thus reverses 
sign, from positive to negative, vanishing at Nf = N 
where 



N 



f,b2z 



2T f (5C A + 3Cf 



(3.5) 
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(The subscript biz stands for the condition that bi is 
zero). Since 



Nf,, 



■N 



/,62a - 



3C a (11C/ + 7Ca) 
4T f {3C f +5C A ) 



> 



(3.6) 



i.e., Nf^ max > Nf t b2z, it follows that there is always a 
nonvacuous interval in the variable Nf where the theory 
is asymptotically free and the two-loop {21) beta function 
has an infrared zero, namely 



Nf, b 2z <Nf < Nf,, 



This zero occurs at 



OtIR.21 



Airbi 



(3.7) 



(3.8) 



and is physical for 62 < 0. Explicitly, for the fundamental 
representation, 



N f,b2z,fv 



UN J 



13A 2 - 3 



and 



Ot-IR,2£,fund 



4n(llN-2Nf 



-34N 2 + Nf(13N - 3N- 1 ) 



(3.9) 



(3.10) 



The sizes of ^-loop contributions arc determined by 
(a/ir) e multiplied by corresponding powers of various 
group invariants. Illustrative values of am t 2t,fund are 
given in Table HVl for N = 2,3,4 and the subset of the 
interval (|3.7|) for which aiR^tjund is not so large as to 
render the two-loop perturbative calculation obviously 
unreliable. Here and below, when a and 7 values are 
listed without an explicit R, it is understood that they re- 
fer to the fundamental representation. Examples of cases 
that we do not include in the table because the two-loop 
result cannot be considered reliable include the follow- 
ing (with formal values of aiR^ijund listed): N = 2, 
Nf = 5, where aiR^tjund = 11-4; N = 3, Nf = 9, 
where chir^i, fund — 5.2; and N — 4, Nf = 11, 12, where 
0!IR,2£Jund = 14, 3.5. 

For reference, the estimate in Eq. (|9.3p of a cr from the 
analysis of the Dyson-Schwinger equation for the fermion 



propagator, in the one-gluon exchange approximation, 
takes the form in Eq. (|9.4p for a fermion in the fun- 
damental representation. This has the respective values 
1.4, 0.79, and 0.56 for N = 2, 3, 4, respectively (where 
we quote the results to two significant figures but do not 
mean to imply that they have such a high degree of accu- 
racy). Setting a cr = otiR t 2l yields the resultant estimates 
of Nf lCr , which, rounded to the nearest integers, are 8, 
12, and 16 for these values of N. We denote these as 
f3DS estimates since they combine a calculation of cxir 
from the perturbative two-loop j3 function with the (one- 
gluon exchange approximation to the) Dyson-Schwinger 
equation. 

As Nf approaches its maximum value, Nf. max , al- 
lowed by the constraint that the theory be asymptot- 
ically free, 62 reaches its most negative value, namely 
b 2 = -C A (7C A + 11(7/). Clearly, for Nf values such 
that 62 is only negative by a small amount and otiR^t 
is large, the perturbative calculation is not reliable. As 
Nf increases further in the range Q3.1[) and olir^I de- 
creases, the calculation becomes more reliable. In Table 
HT1 we list the numerical values of Nf^2z for some illus- 
trative values of N. At the two-loop level, depending on 
whether aiR^i is smaller or larger than a critical value 
for fermion condensation, this is an exact or approxi- 
mate infrared fixed point (IRFP) of the renormalization 
group for the gauge coupling. The existence of such an 
IRFP is of fundamental importance in determining how 
the theory evolves from the ultraviolet to the infrared 
[40| . In particular, as mentioned above, this determines 
whether, as the scale /i decreases sufficiently to a scale 
A (depending on the group G and the fermion content), 
a grows to a large enough size to produce fermion con- 
densates or, on the contrary, the coupling never gets this 
large and the theory evolves into the infrared in a chirally 
symmetric manner, without ever producing such fermion 
condensates. Note that in the former case, the fermions 
involved in the condensates get dynamical masses of or- 
der A and are integrated out of the effective low-energy 
field theory applicable for scales /1 < A, so that the fur- 
ther evolution into the infrared is governed by a different 
beta function. 

It is useful to observe how rapidly the numbers Nf^2z 
approach their large- A values. The number Nf^2z,fund 
has the large- iV expansion 
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N 



34 
13 



102 



306 



(13A) 2 (13) 3 A 4 
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A 6 ) 



N 



2.615- 



0.60355 0.1393 



N 2 



A 4 



O 



N 6 ) 



(3.11) 



As is evident from Table |TT1 the values of Nf^zjund approach the leading asymptotic form for moderate A, as a 
result of the fact that the subleading term in Eq. (|3.11[) is suppressed by 1/A 2 . 
It is of interest to consider the 't Hooft large-A limit, where 



A — > 



with aN fixed. 



(3.12) 



In a theory with fermions in the fundamental representation, in order for them to have a non-negligible effect in this 
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limit, one considers the simultaneous Veneziano limit 

Nj oo with r=^- fixed. (3.13) 

In the combined limit of Eqs. (|3.12[) and (|3.13p . the range of r satisfying the requirement of asymptotic freedom and 
the condition that 62 < so that the two-loop beta function has an IR zero is [41j 

^ < r < y , i.e., 2.615 <r< 5.5. (3.14) 



C. Coefficient bs and Three-Loop Behavior of the Beta Function 



The three-loop beta function coefficient 63 is a quadratic function of Nf with positive coefficients of its N® and N 2 



terms and a negative coefficient of its Nf term. Hence, regarded as a function of the formal real variable Nf, 
positive for large negative and positive Nf, and positive at Nf = 0. The derivative of b 3 with respect to Nf is 



it is 



db 3 
dN f 



Tf 



1415 
~27~ 



C\ 



205 



C A C f 



-2Cj 



T 



§8 r ,316 
jC f + — C A 



(3.15) 



For the fermions representations R that we consider here, for small values of Nf, this derivative db 3 /dNf is negative, 
so that in this region of Nf, 63 decreases from its positive value at Nf — as Nf increases. Because 63 is a quadratic 
polynomial in Nf, the condition that it vanishes gives two formal solutions for Nf, namely 



N 



f,b3z,± 



(U15C A 



615C A C f - 54C? ± Z^Frm 



4Tf(79C A + 66C f ) 



1,2, 



(3.16) 



and 



F Rb3 = 122157C1 + 109578(7^(7/ + 25045C*i<7 : 



7380CUC5 



3240 



(3.17) 



Given that Fr^ > 0, as is the case here, so that the 
values Nf t b3z,j are real, it follows that b 3 is positive in 
the intervals Nf < A/ jb3Zj _ and Nf > Nf^ b3z ^ + and neg- 
ative in the interval Nf^zz,- < Nf < Nf t b3z,+- The 
value Nf t bd,z,+ and the neighborhood of Nf values in the 
vicinity of Nf^ 3z + are not of interest here because they 
are larger than the maximal value Nf^ max allowed by the 
requirement of asymptotic freedom, 



N 



f,b3z,+ 



> N 



f,max 



(3.18) 



Thus, 63 only changes sign once for Nf in the asymptot- 
ically free interval < Nf < Nf_ max . As Nf approaches 
Nf m ax from below, 63 decreases to a negative value given 
by' 



Ca 
24 



1127C A + UC f (UC A -3Cf) . 

(3-19) 

For fermions in the fundamental representation, this is 



(6»)^,^ = -^ + ^. + ^.(3.20) 



As is clear from Table |TT1 for this case 

Nffi3z,l < Nf :b 2z ■ 



(3.21) 



We noted above that any physically reliable zero of the 
beta function must be present already at the level of the 
two-loop beta function, since this is the maximal scheme- 
independent part of this function. Hence, in analyz- 
ing such a zero for the case under consideration where 
the fermions transform according to the fundamental 
representation of SU(AT), we only consider the interval 
(I3.7[) . Combining this fact with our results (13.211) and 
(I3.18p . it follows that b 3 is negative throughout all of 
the interval p.7[) of interest here. For this fundamental- 
representation case, the Nf^zj with j = 1,2 have the 
large- W expansions 



Nf,b3z,l 

and 

Nf,b3z,2 



N 



1.911 



0.3244 0.06844 



N 2 



N 4 



+ 



v 6 ) 



(3.22) 



= N 



13.348 



1.667 0.3978 



N 2 



A 4 



O 



(a 6 ) 



(3.23) 



Here, N ft7nax = 5.5A. 

At three-loop order, the equation for the zeros of the 
beta function, aside from the zero at a = 0, is 61 + 62a + 
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b^a 2 — 0. Formally, this equation has two solutions for a 
and hence for a, namely 



2tt 

ai0z,3£,± = -r~ 
03 



62 ± A/61 - 4&163 



(3.24) 



Since 62 must be negative in order for the beta function 
to have a scheme-independent infrared zero, and since 
for fermions in the fundamental representation we have 
shown that 63 < in the relevant interval Q3.7p , we can 
rewrite this equivalently as 



2tt 



46x16s 



(3.25) 



In order for a given solution to be physical, it must 
be real and positive. As is evident from Eq. (|3.25p . the 
solution corresponding to the the + sign in Eq. (|3.24[) 
(i.e., the — sign in Eq. (|3.25j) ) is negative and hence 
unphysical. Thus, there is a unique physical solution 
for the IR zero of the beta function to three-loop order, 
namely 



OClR,U = apz,3£, 



(3.26) 



Illustrative values for this IR zero of the beta function at 
three-loop order are listed in Table IIV1 

For an arbitrary fermion representation for which /3 
has a two-loop IR zero, we observe that the value of this 
zero decreases when one calculates it to three-loop order, 



(3.27) 



This can be proved as follows. We have 



2tt 
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26 1 |6 3 |+6|-|6 2 K/62 + 46x16. 



(3.28) 

The expression in square brackets is positive if and only 
if 



(26l|6; 



3|+6 2 ) 2 ^ 



b\{b\ +46x|6 3 |) >0 



(3.29) 



But the difference in fl3.29[) is equal to the positive- 
definite quantity b\b\, which proves the inequality (|3.27[) . 
This inequality is evident in Table IIV1 



D. Coefficient 64 and Four-Loop Behavior of j3 



The four-loop beta-function coefficient, 64, was calcu- 
lated in Ref. [29|]. We next analyze its behavior as a 
function of Nf and the result four-loop IR zero of the 
beta function. The coefficient 64 is a cubic polynomial 
in Nf which has positive coefficients of its Nf, and Nf 
terms. Hence, regarded as a function of the formal real 
variable Nf, 64 is negative for large negative Nf, positive 
for Nf — 0, and also positive for large positive Nf. For 
fermions in the fundamental representation, the deriva- 
tive at Nf = is 



/485513 2 ^(3)V 3 +^-^C(3)V+f-^ + ^C(3)V^-?^ 3 (3-30) 



dN f V 1944 9 ' V 1944 9 ^ v ' \ 216 9 



where £(z) is the Riemann zeta function, 

00 

C00=£- (3-31) 

n— 1 

and £(3) = 1.20205690.... This derivative is negative for 
all N. (In the complex N plane, it has six zeros at three 
complex-conjugate pairs of N values.) It follows that, 
again as a function of the formal real variable Nf, 64 
has a local maximum at a negative value of Nf and then 
decreases through positive values as Nf increases toward 
and passes through into the interval of physical values. 

The detailed behavior of 64 in the physical asymp- 
totically free interval < Nf < Nf :inax depends on 
N. In particular, one may determine the value of Nf 
where 64 has a minimum and whether 64 has any ze- 
ros for positive Nf. For SU(2), 64 decreases to a mini- 



mum positive value as Nf ascends through the approxi- 
mate value Nf — 5.8, and then increases monotonically 
for larger Nf, so that it is positive-definite for all non- 
negative Nf, in particular, the asymptotically free region 
< Nf < 11. For SU(3), 64 is also positive-definite 
for all non- negative Nf, reaching a local minimum as 
Nf ascends through a value of approximately 8.2 and 
then increasing monotonically for larger Nf. However, 
for SU(4), 64 is positive for < Nf < 9.51, negative for 
the interval 9.51 < Nf < 11.83, and positive again for 
Nf > 11.83, with zeros at Nf ~ 9.51 and Nf ~ 11.83. 
We list these zeros of 64 as a function of Nf in Table [IT] 
(Again, we recall that the physical values of Nf are, of 
course, restricted to non- negative integers.) Thus, this 
reversal of sign occurs in the interval of interest here, 
< Nf < 22, where the SU(4) theory is asymptotically 
free. For SU(5), 64 behaves in a manner qualitatively sim- 
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ilar to the SU(4) case; it is positive for < Nf < 11.18, 
negative in the interval 11.18 < Nf < 15.18, and posi- 
tive for larger values of Nf, vanishing at Nf ~ 11.18 and 
Nf ~ 15.18. Thus, again, 64 reverses sign in the region 
< Nf < 22.5 where the SU(5) theory is asymptotically 
free. Thus, in contrast with 62 and 63, which are negative 
throughout the interval of Nf of interest (and 61 , which 
is positive), 64 can, for N > 4, vanish and reverse sign in 
this interval. 

At the four-loop level, the equation for the zeros of the 
beta function, aside from a — 0, is the cubic equation 



bi + b 2 a + 63a + 64a' 







(3.32) 



This equation has three solutions for a and hence for 
a, which will be denoted ap Zj u,j, j = 1,2,3. Since the 
coefficients be are real, there are two generic possibilities 
for these three roots, namely that they are all real, or that 
one is real and the other two form a complex-conjugate 
pair. The properties of the roots are further restricted 
by the asymptotic freedom condition that b\ > 0, the 
existence of a two-loop IR zero, which requires that 62 < 
0, and the fact that, as we have shown, for the relevant 
range (|3.7p of Nf, where these conditions are met, 63 < 0. 
As is evident in Table IIV1 we find that for the values of 
N and Nf that we consider, the roots of Eq. (|3.32l) arc 
real. For all of the values of N and Nf where there is a 
reliable two-loop value for an IR zero of the beta function 
(i.e., where it does not occur at such a large value of a 
as to render the perturbative calculation untrustworthy) , 
one of these roots is negative and hence not physical, one 
of them, namely the minimal positive one, is the physical 
IR zero, which we will denote am 44 — am,4t/ '(4tt), and 
there is a third root at a larger positive value. This third 
root, denoted au,u = <^ai,u/ (47r), is not relevant for our 
analysis, since the initial value of a at a high energy scale 
/1 is assumed to be close to zero, so that as the scale n 
decreases, a increases and approaches the (positive) zero 
of the beta function closest to the origin, namely am u 
42]. 

It is straightforward to display the analytic expressions 
for the root am,u, but we shall not need this for our 
analysis. We list numerical values for am,Ai for various 
values of N and Nf in Table IIVI For completeness, we 



note the specific sets (N,Nf) where am,2( is so large 
that we consider the analysis via the perturbative beta 
function unreliable: these are (N,Nf) — (2,6), (3,9), 
(4,11), and (4,12). 



E. Estimates of Zeros of the Four-Loop Beta 
Function via Pade Approximants 

For the beta function, or more conveniently, 
the reduced function with the pref actor removed, 
Sj=o 3!_1 bjaP -1 , it is useful to calculate and analyze 
Pade approximants, since these provide closed-form ex- 
pressions that, by construction, agree with the series to 
the maximal order to which it is calculated. The ex- 
pansion for j3 a to t = 4 loop order can be used in two 
ways. First, one can simply solve the cubic equation 
Pa = bi + b2a+b 3 a 2 +b4d 3 = and obtain the three roots, 
one of which is the root of interest, giving the IR zero. 
Secondly, one can calculate Pade approximants, e.g., the 
[2,1] and [1,2] approximants, and determine their zeros. 
The [1,2] Pade approximant has a single zero at 



a /3z,4£,[l,2] 



a iRAi,[i,2] _ h(bib 3 - bl) 



4tt b 3 2 ~2b 1 b 2 b 3 + b 2 1 b i ) 



(3.33) 



Taking into account the fact that 62 and 63 are negative 
in the relevant interval (I3.7[) . this can be rewritten as 



a Pz,Al,[l,2\ 



a /i?,4£,[l,2] 
47T 



6i(6i|6 3 |+&i) 



|62| 3 + 26 1 |6 2 ||o 3 |-o 1 &4 



The two zeros from the [2,1] approximant are 



(3.34) 



3z,4£,[2,l],± 



6263 - 6164 ± 



(b2b 3 ~hb 4 ) 2 -4bMH- « 



1/2 



2(6204 - b 2 3 ) 



(3.35) 



Taking account of the fact that 62 and 63 are negative in 
the relevant interval (|3.7I) . this can be rewritten as 



a /3z,4£,[2,l],± 



6164 - MN T (MN - 0164) 2 + 4641631(61 + I62I64) 



1/2 



2(|6 2 |6 4 
I 



(3.36) 



The expression in Eq. (|3.36p with the — sign in front 
of the square root is negative and unphysical, while the 
expression with the + sign in front of the square root 
yields the estimate of the IR fixed point, as en/ R,je, \2,i] = 
4:7ra/3z.4£,[2,i]- As is evident from Eqs. (|3.33[) and (|3.35l) . 



the zeros of the [1,2] and [2,1] Pade approximants incor- 
porate information on j3 up to four loops. One readily 
verifies that in the limit 64 — > 0, the zero of the [1,2] Pade 
reduces to the two-loop result a = —61/62, and the two 
zeros of the [2, 1] Pade reduce to those obtained from the 
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three-loop beta function, (|3.26[) . We list the values of 
a IR obtained from the zeros of the [1, 2] and [2, 1] Pade 
approximants to the four-loop beta function for the case 
of fermions in the fundamental representation in Table 

El 

From our calculations of am at the three- and four- 
loop level for SU(A) with fermions in the fundamen- 
tal representation, we can make several remarks. Al- 
though n-loop calculations of the beta function for n > 3 
loops are scheme-dependent, the results obtained with 
the present A/5 scheme provide a quantitative measure 
of the accuracy of the scheme-independent two-loop re- 
sult. For a given A, as Nf increases above the mini- 
mal value Nf t b2zt where the IR zero first appears, and as 
the resultant aiR^e decreases to values < 1, the differ- 
ence between oljr^i and the higher-loop values aiR. n i for 
n = 3, 4 decrease. As is evident from Table [IV] the value 
of ociR,nl generically decreases as one goes from n — 2 to 
n — 3 loops and then increases by a smaller amount as 
one goes from n — 3 to n = 4 loops, so that aiR^e is 
smaller than aiR^e- In the same region of Nf values 
such that am,2i is reasonably small, the values obtained 
via the [1,2] and [2,1] Pade approximants to the four-loop 
beta function are close to those obtained from the zeros 
of this beta function itself. 



IV. EVALUATION OF THE ANOMALOUS 
DIMENSION 7m AT THE INFRARED ZERO OF /3 

In this section we evaluate the anomalous dimension of 
7 = 7 m , calculated to the n-loop order in perturbation 



theory, at the (approximate or exact) IR zero of the beta 
function to this order, ajR, n ^, for n = 2, 3, 4. We denote 
these as j n e(oiiR t ne)- We focus here on general results 
and their application to the case of fermions in the fun- 
damental representation, and discuss higher-dimensional 
representations in subsequent sections. In general, this 
anomalous dimension must be positive to avoid unphys- 
ical singularities in fermion correlation functions. 

A running fermion mass, S(fc), that is dynamically 
generated at a scale A, decays with Euclidean momen- 
tum k > A like 

£(fc)~A(-J (4.1) 

up to logs. Since for k > A, the running coupling a 
is smaller than the critical value a.R, cr and there is no 
spontaneous chiral symmetry breaking, it follows that 
S(fe) must decrease toward zero as fc/A — > oo. In turn, 
this implies that j m < 2. Hence, a physical value of 7 m 
must lie in the range 

< 7m < 2 . (4.2) 

For values of Nf such that the theory evolves into the 
infrared in a chirally symmetric manner, so that the IR 
zero of the beta function is exact, the same upper bound 
follows from a related unitarity consideration [43j. 

Using the two-loop result for 7 and evaluating it at 
the two-loop value of the IR zero of the beta function, 
we have 



Cf(llC A - AT f Nf)(455C 2 A + 99C A C f + (180C/ - 248C A )T f N f + 80T 2 AT 2 ) 
12(-17C 2 +(10CU + &Cf)T f NfY 



For the fundamental representation, this is 

(N 2 - 1)(1L/V - 2A^)(1009A 3 - 997V - (158A 2 + 90)A^ + AONNj) 



l2l{aiR,2l) 



12(-34N 3 + (137V 2 - 3)AT/) 2 
I 



(4.3) 



(4.4) 



We list numerical values of 7(0^,2*) in Table IVT1 for the 
illustrative values N = 2, 3, 4 and, for each N, a set of 
Nf values in the range (I3.7[) . For sufficiently small Nf > 
Af/,b2z in each case, aiR^e is so large that the formal 
value of ^tip-iR^i) is larger than 2 and hence unphysical; 
we enclose these values in parentheses to indicate that 
they are unphysical artifacts of a perturbative calculation 
at an exessively large value of a. 



In the large- A, large- N f limit of Eqs. (pTL2l and (|3TT3"|) 



with r ee Nf/N, Eq. (|4.4p reduces to 



j2e(aiR^e) = 



(11 - 2r)(1009 - 158r + 40r 2 ) 



12(-34+ 13r) 2 



-01 



A 2 ) ' 
' (4-5) 

For r — 4 corresponding to the asymptotic value of 
N LcrJ und in Eq. (HU), j2i(a IR} 2e) = 113/144 ~ 0.785, 
which is the same as the large- A limit of Eq. (14.71) . 

One may evaluate ^lio-iR^i ) at Nf equal to the value 
Nf^cr.fund predicted by the one-gluon exchange (ladder) 
approximation to the Dyson-Schwinger equation for the 
fermion propagator, given in Eq. (|9.8p . This is somewhat 
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formal, since these values of Nf tCr j un d are not, in general, 
integers and hence not actually physical; for example, 
N f , crJund = 7.86, 11.91, 15.94 for N = 2, 3, 4). This 
procedure yields the result 

565A 4 - 7067V 2 + 225 
72e{ai R , 2 i,Nf, cr Jund) - 144 (JV2- i)(5JV 2 -3) ' ( ' 

For the illustrative cases N — 2, 3, 4, this anomalous 
dimension takes the values 0.88, 0.82, and 0.80, respec- 
tively. As N — > oo, Eq. (|4.6p has the expansion 



, . 113 11 



144 ' 40 A 2 ' \N 4 



1 



(4.7) 



Since the estimate (|9.8p is close to 4A even for the small- 
est value, N = 2, and asymptotically approaches 47V as 
N — > oo, it is worthwhile to compare the above values 
of 7, viz., 0.88, 0.82, and 0.80 for N = 2, 3, 4, with 
l2i{oiiii.2t) evaluated at the nearest physical, integer val- 
ues of Nf, namely N f = 8, 12, 16 for N = 2, 3, 4. 
This procedure yields j2i(c(iR,2e) = 0.75, 0.77, 0.78, as 
recorded in Table [VTl To within the strong-coupling the- 
oretical uncertainties of these calculations, these values 
are mutually consistent. 

A closely related approach is to evaluate the two-loop 
expression for j m at a = a cr> R, where a crt u is the es- 
timate of the critical coupling for fermion condensation 
obtained from the one-gluon exchange approximation to 
the Dyson-Schwinger equation, given in Eq. (|9.3[) . If 
one were to use the first term in Eq. (|2.5[) by itself, one 
would obtain, to this order, j m — ci{a cr ^/i:) — 1/2, for 
an arbitrary fermion representation R. The second term, 
c 2 {oi C r,R/^) 2 , does depend on R. Combining these, i.e., 
evaluating ci(a/ir) + ^(a/ir) 2 at a = a cr ,Rt we find, for 
a general fermion representation R, 



l2i{a. 



cr,R) 



97C A + 225C> - 201/ TV/ 
432C7 



(4.8) 



For consistency, one would then substitute Nf equal 
to the estimate Nf iCrt R from the same Dyson-Schwinger 
analysis, given in Eq. (|9.7I) . For a general representation 
this yields 



2\C\ + 128C A C f + 225(7f 
144(7/ (CU + 3(7/ 



l2l(ctcr,R\ Nf = Nf t0rtR ) 

1 -t-H /'ll \ - ( ! 

For the fundamental representation, this reduces to the 
same result as was obtained in Eq. (I4.6|) . We have also 
evaluated the three-loop result for 7 at the three-loop 
value of the IR zero of the beta function, which we denote 
as 73£(a/fl.3f), and the four-loop result for 7 at the four- 
loop value of the IRFP, which we denote as "1ai{o.ir,u)- 
We list the resultant values in Table PvTl 

From our calculations of j m for the case of fermions 
in the fundamental representation, we can make sev- 
eral observations. Although computations of am^i and 
lni{ctiR,ni) are scheme-dependent for n > 3 loops, they 



provide a useful measure of the accuracy of the lowest- 
order results. As was the case with the position of am tn £ 
itself, we find that, for a given N and for Nf reasonably 
well above Nf^2z so that the perturbative calculation of 
otiR t ni is not too large, the value of ^nuipHR^t) generi- 
cally decreases as one goes from n = 2 to n = 3 loops. 
Some of this decrease can be ascribed to the decrease in 
aiR,n£ going from (n = 2)-loop to (n = 3)-loop order. At 
the four-loop level, ^u{ctiR,u) tends to be smaller than 
"tu{&iR,u) f° r values of Nf from Nf^z to values of Nf 
slightly above the middle of the range (|3.7|) . while for 
values of Nf in the upper end of this range, jui&iR^i) 
is slightly larger than 73^(077?, 3^). In general, for the 
values of Nf where am is sufficiently small that the cal- 
culation may be trustworthy, the value of the anomalous 
dimension evaluated at the IR zero of the beta function 
(both calculated to n-loop order) ~i n i{a-iR,ni)i 1S some- 
what smaller than unity. 

Several recent high-statistics lattice simulations have 
been carried out on an SU(3) gauge theory with a vary- 
ing number Nf of fermions in the fundamental represen- 
tation in the range 6 < N f < 12 [ll|-[il, 0, [H|. This 
work has yielded evidence for a regime of slowly running 
gauge couplings for Nf < 12, consistent with the pres- 
ence of an IR zero of the beta function, in agreement 
with the earlier continuum estimates in Ref. 1]. Ref. 
[Ill ] also found a considerable enhancement of ('0'0)//|> 
in the SU(3) theory with Nf = 6. Further lattice sim- 
ulations and analysis of data should yield values of 7 m 
that can be compared with our higher-loop calculations 
in this paper. A preliminary study of the SU(2) theory 
with Nf = 6 fermions has also been reported [211 ]. 



V. ADJOINT REPRESENTATION 

In this section we analyze the SU(iV) theory with Nt 
copies of a Dirac fermion, or equivalently, 2 AT/ copies of 
a Majorana fermion, in the adjoint representation. For 
this case, the general expression for the maximal value 
of Nf allowed by the requirement of asymptotic freedom, 
Eq. ([321), reduces to 



N 



f ,max ,ad 



11 
T 



(5.1) 



i.e., restricting Nf to the integers, Nf^ max = 2. The gen- 
eral expression in Eq. (j3.5[) for the value of Nf at which 
62 changes sign from positive to negative with increasing 
Nf reduces to 



N ftb2Ztadj = ^ = 1.0625 



(5.2) 



Hence there is only one (integer) value of Nf, namely 
2 Dirac fermions (equivalently, Nf = 4 Majorana 



Nf 

fermions), for which the theory is asymptotically free and 
has an IR zero of the two-loop beta function. This zero 
occurs at 

27r 1 257 

0tIR,2l,adj = 5JY - ~TT for N f = 2 ■ ( 5 - 3 ) 



N 
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Specializing the general formula for the critical cou- 
pling a CTi R from the one-gluon exchange approximation 
to the Dyson- Schwinger equation, Eq. (|9.3j) (see Ap- 
pendix II) for the present case where R is the adjoint 
representation, one obtains a cr ^ a dj = n/(3N). Formally 
setting a>iR,2£,adj = a> a dj,cr yields the corresponding esti- 
mate for the critical number iV/ )Cr = 83/40 = 2.075. This 
may be rounded off to the nearest integer, giving Nf jCr for 
the adjoint representation. In view of the theoretical un- 
certainty in such an estimate, due to the strong-coupling 
nature of the physics involved, an SU(iV) gauge theory 
with Nf = 2 adjoint fermions could be either slightly in- 
side the chirally broken, confined side of Nf >cr or slightly 
on the other side, where the theory is chirally symmet- 
ric and the evolution into the infrared is governed by an 
exact conformal IR fixed point. 

For the present case of Nf = 2 fermions in the adjoint 
representation of SU(JV), the coefficients of the beta func- 
tion are h = N,b 2 = -107V 2 , b 3 = -(101/2)iV 3 , and 



5977 



54 
176 
"IT ~ 



N 4 
176 



+ 27^ 



448 
~27 



N 



128 



N 



C(3) 

416 
~3~ 



-N 4 



N 



3A 2 
124 
~9 

128 
~3~ 



208 
~~9~ 

128 

' TV 2 



JV' 



3232 



N 



(5.4) 



(The term in C3 proportional to £(3) and the terms in 
c 4 proportional to £(4) and £(5) vanish for the adjoint 
representation for arbitrary Nf.) 

Evaluating the two-loop expression in Eq. (|4.3p for 7 m 
at the IR zero of the beta function, also calculated at the 
two- loop level, aiR^i^dj, we obtain 



(11 - 4iV / )(277 - MN f + 40JV 2 ) 



6(-17+ 167V/) 2 



so that for the Nf = 2 case of interest here, 



(5.7) 



41 



I2i,ada{am at ,ada) = ^ =0.820 for Nf = 2 . (5.8) 

It is also of interest to evaluate the two-loop 7 m at 
the value of a cr from the one-gluon exchange (ladder) 
approximation to the Dyson- Schwinger equation. Spe- 
cializing the general result in Eq. (|4.8p to the adjoint 
representation with Nf = 2 yields 



47 

72^ (oicr,adj) = 72 ~ 0-653 



(5.9) 



Evaluating the three-loop result for 7 m at the IR zero 
of the beta function calculated at the three-loop level, 
o-iR,3t,adj, for the Nf = 2 case of interest, we obtain 



jM,adj(®iR,M,<*dj) = 0.543 for N f = 2 



(5.10) 



At the three-loop level, the beta function has two zeros 
away from the origin, one of which is negative and un- 
physical, and the other of which is the physical IR zero, 



1017V 



0.9180 
~N~ 



(5.5) 



At the four-loop level, the beta function has three zeros 
away from the origin, one of which is negative and un- 
physical, one of which is the four-loop IR zero, otiR t u,adg, 
and the third of which, denoted au.u, occurs at a larger 
positive value of a and is not relevant to our study, since 
it is not reached by evolution of the coupling starting at 
small a for large \i. We list the numerical values of these 
zeros in Table Ell 

The coefficients q in Eq. (|2.5[) for 7 for this case are 
ci = 3N/(2n), c 2 = (lliV 2 )/(87r 2 ), and c 3 = -iV 3 /(27r 3 ), 
with C4 given by 



4- 
7T C4 



1817 
512 



-N 



1 N 3 _ JLjv 2 

96 192 



7 _ _ 1 1 

— iVH 1 

96 32 4iV 



C(3) 



155 , 5 o 
■N N 3 



128 



64 
15 



16iV 3 
25 
128' 



-N 



45 



35^_ 15 

64 ~ 64 _ 87V + 32iV 3 



(5.6) 



which is again independent of N. At the four-loop level, 
the value of j4e !a dj(<xiR,4£,adj) does depend slightly on N. 
We list the values of these anomalous dimensions in Ta- 
ble IVIIII The most recent simulations of a lattice gauge 
theory with SU(2) gauge group and Nf — 2 fermions in 
the adjoint representation report j m = 0.49 ± 0.13 [23|. 
This is in agreement with the calculations of j m here at 
the three- and four-loop level, to within the uncertainties 
of the respective calculations. 



VI. SYMMETRIC AND ANTISYMMETRIC 
RANK- 2 TENSOR REPRESENTATIONS 

In this section we consider the SU(iV) theory with Nf 
fermions in the symmetric or antisymmetric rank-2 rep- 
resentation, denoted S2 and A2. Since a number of for- 
mulas are similar for these two cases, we will often give 
these in a unified way for both cases, denoted T2 (for 
rank-2 tensor representation), with ± signs distinguish- 
ing them. For S2, our analysis applies for any N, while 
for A2, we restrict to N > 4, since the A2 representation 
is the singlet for SU(2) and is equivalent to the conju- 
gate fundamental representation for SU(3). Note that 
for SU(4), the A2 representation is self-conjugate. Also, 
since for SU(2) the S2 representation is the same as the 
adjoint representation, which has already been analyzed, 
we only consider the illustrative values N = 3, 4. 
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For the two T2 cases, the general expression for the 
maximal value of Nf allowed by the requirement of 
asymptotic freedom, Eq. (|3.2|) . reduces to 



N 



UN 



f,max,T2 



2(N±2) 



(6.1) 



where the ± refers to S2 and A2, respectively. As N 
increases from 2 to oo, Nf, maXj s2 increases monotonically 
from 2.75 to 11/2 = 4.5, and as N increases from 3 to oo, 
Nf, m ax,A2 decreases monotonically from 16.5 to the same 
limit, 4.5. The physical values of Nf lTnax in both cases 
are the greatest integral parts of these rational numbers. 

For these representations, the general expression in Eq. 
p.5[) for the value of Nf at which the beta function co- 
efficient 62 changes sign from positive to negative with 
increasing Nf takes the form 



N 



17 N 2 



f,b2z,T2 



(N ±2)(8N ±3-6iV-i) 



As a consequence of the general inequality (|3.6I) , it follows 
that N f , b2z ,S2 < N f and Nf,b2 Z ,A2 < Nf ,max,A2 ■ 

For N — 2, the S2 representation is just the adjoint rep- 
resentation, so we only consider the illustrative values 
N = 3, 4. The respective intervals Nf ! b2z,S2 < Nf < 
Nf, m ax,S2 for which the SU(AT) gauge theory is asymptot- 
ically free and has an IR zero of f5 are 1.06 < Nf < 2.75 
for N = 3 and 1.22 < Nf < 3.30 for N = 4. These 
ranges imply that the only physical integral values of Nf 
satisfying these conditions are Nf = 2, 3 for both SU(3) 
and SU(4). 

For large N, iV/ 62«,T2 nas the series expansion 



N 



17 323 6137 103547 



f,b2z,T2 



2 3=F 2 6 AT 



2 9 N 2 * 2 12 N 3 



-0 



I 



As N increases from 2 to oo, A//,62«,S2 increases mono- 
tonically from 17/16 = 1.0625 to 17/8 = 2.125, and as 
N increases from 3 to 00, Nf t b2z,A2 decreases monoton- 
ically from 8.05 to the same limit, 2.125. This limit is 
twice the (A^-independent) value of Nf^2z,adj — 17/16 
for the adjoint representation. Thus, for large N, the 
range (|3.7|) where the SU(iV) theory with Nf fermions in 
the S2 or A2 representation is asymptotically the same 
for both, namely, 17/8 < Nf < 11/2; restricting Nf to 
physical, integer values, this range consists of the three 
values Nf = 3, 4, 5 

For our further discussion we assume that Nf is in the 
range Nf t wsz,T2 < N f < Nf^ max , T2 where the theory is 
asymptotically free and the two-loop beta function has 
an IR zero, for the respective cases S2 and A2. This zero 
occurs at the value 



2tt(UN -2N f {N ±2) 



-17iV 2 + N f (8N 2 ± 197V =F 12AT" 1 ) 



(6.4) 



Evaluating Eq. (|9.3[) from the ladder approximation to 
the Dyson- Schwinger equation for the fermion propaga- 
tor for S2 or A2 representations yields the estimate (|9.6j) 
given in Appendix II. The resultant (3DS estimates for 
Nf_ cr in the case of the S2 representation and TV = 2, 3, 4 
are Nf >crj s2 — 2.1, 2.5, 2.8, respectively. For the A2 rep- 
resentation with N = 4, one has -/V/ cr ,A2 = 8.1. 

The two-loop expression for the anomalous dimension, 
evaluated at a = aiR2e.,T2, is 



l2l,T2(aiR.2e,T2) — 

(N ± 2)(N T 1) [llAT - 2(N ± 2)N f 



N{5MN 2 ± 99N - 198) + (-34N 3 ± 22A 2 =f 360) + 20N (N ± 2) 2 Nf 



12 



17N 3 + {N± 2)(8N 2 ±3N- 6)N f 



(6.5) 



r 



We list values of "i2t,S2 {&ir,21, S2) for N = 2, 3, 4 in Table 
Poland values of 72^2 {ohr,21,A2) for N = 4 in Table |XTI1 
It is also of interest to evaluate the two-loop expression 
for 7 at the estimated a = a cr ,T2- This yields 



72£,T2(a C r,T2) 



322N 2 ± 225iV - 450 - 10N(N ± 2)N f 



432(A^±2)(iV=F 1) 



We list these values in Tables IXII and IXIII 



(6.6) 



Evaluating the two-loop anomalous dimensions at the 
two- loop IR zero of the beta function, ^2i,T2 {oiin,2t,T2), 
for Nf equal to the respective BDS-estimated critical val- 
ues in Eq. (|9.10[) . we obtain (again with T2 and the ± 
signs referring respectively to S2 and A2) 



J21,T2 {<XIR,2£,T2) \Nf—N fiCriT2 = 

374N 4 ± 578AT 3 - 931N 2 qp 900N 



900 



144(AT ± 2) (AT =p 1)(4A^ 2 ± 3AT - 6) 



(6.7) 
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This has the large- TV expansion 

"/2e,T2(aiB„2e.T2)\N f =N JtCr ^ T2 = 

187 17 1 \ 

m T 128iV + °M- ^ 

The leading term has the value 187/288 ~ 0.649. 

From a lattice study of SU(3) gauge theory with Nt = 
2 fermions in the S2 (sextet) representation, Ref. [18[ 
found that this theory is characterized by slow running 
behavior consistent with an (exact or approximate) IR 
fixed point, and further reported that 7 m < 0.6 where 
it was measured. For SU(3), the estimate of a cr ,S2 in 
Eq. (I9.3[) gives a cr ,s2 — it/ 10 = 0.31. Our results for 
the IR zero of f3 and the value of 7 m at this zero for 
N = 3 and N f = 2 are listed in Tables |IX] and EH We 
find that aiji,ni,S2 is approximately 0.84 at n = 2 loop 
level, decreases somewhat to 0.50 at three-loop level, and 
then increases slightly to 0.52 at four-loop level. The 
two-loop result for 7 m is unphysically large, while the 
three- and four-loop values of 7 m at the corresponding 
three- and four- loop IR zeros of (3 are about 1.3 and 
1.4. These are somewhat larger than the values reported 
in Ref. [lH, although in assessing this comparison, one 
must take account of the significant strong-coupling un- 
certainties in our calculation stemming from the fact that 
&ir,S2 ~ 0(1). Our evaluation of the two-loop expres- 
sion for 7 m at the ladder-Dyson-Schwinger estimate of 
a cr ,S2, is 0.65. 

VII. EFFECTS OF NONZERO FERMION 
MASSES 

The global chiral symmetry that is operative if the 
fermions are massless, and the way that it is broken by 
fermion condensates, is well-known, and we do not review 
it here. However, it is worthwhile to comment on the sit- 
uation in which some fermion masses are nonzero. In this 
paper we generally assume that the fermions have zero in- 
trinsic masses in the Lagrangian describing the high-scale 
physics, and the only masses that they acquire arise dy- 
namically if they are involved in condensates that form as 
the gauge interaction becomes sufficiently strongly cou- 
pled in the infrared. This is a well-motivated assumption 
if the vectorial gauge theory arises as a low-energy effec- 
tive field theory from an ultraviolet completion which is 
a chiral gauge theory. In turn, this is natural if the latter 
theory becomes strongly coupled, since it can then form 
fermion condensates that self-break it down to the vecto- 
rial subgroup symmetry. However, one may also choose 
to focus on the vectorial gauge theory as an ultraviolet- 
complete theory in itself. In a vectorial gauge theory, an 
intrinsic (bare) mass term for a fermion ip, C rn — — m^>, 
is allowed by the gauge invariance. Hence, one may con- 
sider a more general situation in which the fermions may 
have such intrinsic (hard) masses in the high-scale La- 
grangian [44j]. In this case, as the reference scale \i de- 



creases below the value of the hard mass of some fermion 
m/, the beta function changes from one that includes this 
to one that excludes this fermion. If the hard fermion 
masses are small compared with the scale A in the situ- 
ation where the theory confines and breaks chiral sym- 
metry spontaneously, then these hard masses have only a 
small effect. However, if some of the hard fermion masses 
are sufficiently large, then as /i decreases below their scale 
and the corresponding fermions are integrated out of the 
low-energy theory below this scale, this can significantly 
change the infrared properties of the resultant theory. 

In applications of slowly running gauge theories 
to technicolor theories, at the scale Arc where the 
SXJ (Ntc) gauge coupling grows to 0(1) and is influenced 
by the presence of an approximate IR zero of the TC 
beta function, there can also be non-negligible effects 
due to four-fermion operators arising from the higher- 
lying extended technicolor dynamics [To|-j45l. |46|. 
and these can affect the scaling properties of -tptp. Similar 
comments apply for topcolor-assisted technicolor 0, E3] ■ 



VIII. CONCLUSIONS 

In this paper we have studied the evolution of an 
asymptotically free vectorial SU(iV) gauge theory from 
high scales to the infrared taking account of higher-loop 
corrections to the beta function and the anomalous di- 
mension 7 m for fermions in the fundamental, adjoint, and 
rank-2 symmetric and antisymmetric representations S2 
and A2. We have compared our results with lower-order 
calculations. We have shown that, for fixed N and Nf, 
in the range for which the two-loop beta function has an 
IR zero, the value of this zero decreases as one goes from 
the two-loop to the three-loop calculations, and we have 
determined this decrease quantitatively. Going further, 
we have shown that there is a smaller fractional increase 
in the value of this IR zero when calculated to four-loop 
accuracy, with the final four-loop result still smaller than 
the two-loop value. We have analyzed instanton effects 
and have demonstrated that they tend to increase the 
value of the IR zero of the beta function somewhat. A 
major part of our work has been the evaluation of the 
anomalous dimension r y m of ipif) at the IR zero of the beta 
function at the £ = 2, 3, 4 loop levels. This zero is approx- 
imate or exact, depending on whether for a given N, the 
value of Nf is below or above the critical value Nf^ cr be- 
low which there is spontaneous chiral symmetry breaking 
associated with the formation of a fermion condensate. 
We have found that this 7 m at the (approximate or ex- 
act) IR zero of the beta function decreases as one goes 
from two-loop to three-loop order, and that the four-loop 
values also tend to be somewhat less than those at the 
two-loop level. The values that we have calculated for 7 m 
at the IR zero of the beta function tend to be somewhat 
smaller than unity. We have compared our higher-loop 
calculations with results from recent lattice simulations 
and have found general agreement. We believe that the 
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higher-loop calculations reported here should provide a IX. APPENDIX I 

useful reference for comparison with ongoing and future 
lattice measurements. 
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2857 
~54~ 



C\ + +T f N f 



9r 2 205 1415 2 



(TfN f ) 2 



44 158 



(9.1) 



The four- loop coefficient is given in Ref. 29] and is a cubic polynomial in Nf. We note that the coefficients of the 
N® (which is independent of the fermion representation) is positive, and the coefficient of the term is positive for 
an arbitrary fermion representation. 

Our normalizations for the quadratic Casimir and trace invariants of a Lie group are standard. The quadratic 
Casimir invariant Ca(-R) for the representation R is given by Ea=i X^St [DR(T a )]ij[DR(T a )]jk = C2{R)8ik, 
where a, b are group indices, o(G) is the order of the group, T a are the generators of the associated Lie alge- 
bra, and Dn(T a ) is the matrix form of the T a in the representation R. The trace invariant T(R) is defined by 

=i fi) [DR(Ta)]ij [D R {T b )] 3l = T(R)S ab . 

From the calculations of the coefficients of the perturbative expansion of the anomalous dimension j m in the MS 
scheme to four- loop order in Ref. [3(j, we record the three-loop coefficient 



c 3 = 2C 



129 2 



129 

—c f c A 



11413 
108 



Cl + C / T / 7V / (-46 + 48C(3)) 



48C(3)) - ^TfNj 



(9.2) 



We have used the four-loop coefficient C4 from Ref. [3C 



for our calculations, but it is too lengthy to reproduce here. 
I 



A. Appendix II: Beta-Dyson-Schwinger Estimate 

of N f>cr 

In this appendix we collect some useful known formulas 
from the one-gluon exchange (also called ladder) approxi- 
mation to the Dyson- Schwinger equation for the fermion 
propagator 0, 0, IHJ. With this approximation, a so- 
lution of the Dyson-Schwinger in the case of an initial 
zero input mass of a fermion in the representation R of 
the gauge group, yields dynamically generated mass if 
the coupling a(/i) exceeds a critical value a cr r given by 

aaS 



OL C r,R 



3C f 



(9.3) 



For the representations of interes here, this has the form 

27riV 



&cr,fund 



3(N 2 - 1) 



Q-cr,adj 



3N 



TTiV 



Olcr,T2 



3(7V±2)(7V=Fl) 



(9-4) 



(9.5) 



(9.6) 



where, as in the text, the ± signs apply for the S2 and 
A2 representations, respectively. 

For the illustrative values N = 2, 3, 4, this has the 
following values: (i) 1.40, 0.79, 0.56 for the fundamental 
representation; (ii) Setting this equal to the two-loop ex- 
pression for the IR zero of j3 then yields an estimate for 
Nf^ cr to this order, namely 



c A {mc f + nc A ) 

lOT f (C A + 3C f ) 



(9.7) 



For the various representations of interest here, this has 
the explicit forms 



N 



f,cr,fund 



N 



ficr,adj 



2N{50N 2 ~ 33) 
5(57V 2 - 3) 



g = 2.075 
40 



and 



N 



f,cr,T2 



N(83N 2 ± 66iV - 132) 
5(7V±2)(4iV 2 ±3iV-6) 



(9.8) 



(9.9) 



(9.10) 



where, as in the text, T2 refers to S2 or A2 and ± cor- 
responds to S2 and A2, respectively. In the same ladder 
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approximation, one finds 7 m = 1 at a = a cr ^R [lj (which 
also holds for the DS analysis at a UV-stable fixed point 
48]). For the gauge group SU(iV) with the illustrative 
values of N used for the tables, namely N = 2, 3, 4, 
Nf tCr ,fund is equal to 7.9, 11.9, and 15.9, respectively, 
with the large- N form Nf ~ AN. For S2, the symmetric 
rank-2 tensor representation, N = 2, 3, 4, Nf iCr< s2 is 
equal to 2.075, 2.5, and 2.9, increasing toward the limit 
11/2 = 5.5 in the large- N limit. In the case of A2, the an- 
tisymmetric rank-2 tensor repreprescntation, for N = 3, 
the result is the same as for the fundamental representa- 
tion, while for N — 4, one has Nf^ rity A2 — 8.1, and as 
N — > oo, Nf^ cr it t A2 decreases toward the limit 11/2. 

One understands that, a priori, there could be sig- 
nificant uncertainty in these estimates because of the 
strong-coupling nature of the physics involved and the 
one-gluon approximation used for the solution of the 
Dyson- Schwinger equation. Moreover, the DS equation 
analysis is semi-perturbative in the sense that it contains 
polynomial dependence on a, and it neglects nonpertur- 
bative effects associated with confinement and instan- 
tons. However, corrections to the one-gluon exchange 
approximation have been analyzed and found not to be 
too large [2|. Recent lattice simulations for SU(3) are 
in broad agreement, to within the uncertainties, with the 
above prediction of N f . cr ~ 12 [ll|-[ll], [Hj]. Some of the 
success of the f3DS prediction for Nf^ cr may arise from 
the fact that two major physical effects that it ignores 
would shift Nf iCr in opposite directions and hence tend 
to cancel each other out [34J . The DS equation is an inte- 
gral equation, and the standard analysis of this equation 
involves an integration over Euclidean loop momentum 
k from k = to k = oo. If the theory confines, then the 
lower bound for the Euclidean loop momentum should ac- 
tually not be k — 0, but instead k = k m i n . ~ r^ 1 where 
r c is the spatial confinement scale [HJ . The use of k = 
thus overestimates the tendency toward S%SB. Instan- 



tons enhance S%SB [36H39I ] , and the neglect of instanton 
effects leads to an underestimate of the tendency toward 
SxSB. Since these two neglected aspects of the physics 
- confinement and instantons - produce errors that are 
of opposite sign as regards the tendency for S%SB, it is 
plausible that these errors tend to cancel out, so that the 
/3DS method may be reasonably accurate in predicting 

N f,cr- 



X. APPENDIX III - PADE RESULTS 

In this appendix we collect some relevant results on 
Pade approximants. Given a Taylor (or asymptotic) se- 
ries expansion around z = for the function f(z), 



/(*) 



E 



fnZ r ' 



(10.1) 



one can construct a set of \p, q] Pade approximants, 
namely rational functions comprised of a numerator poly- 
nomial of degree p and a denominator polynomial of 
degree q, such that p + q = n max — 1, of the form 
=oPj zj )/{YX=o 1kZ k ). Without loss of generality, 
one can divide numerator and denominator by qo, so that, 
after redefinition of the coefficients, one has 



E 



qkz n 



(10.2) 



The p + q + 1 coefficients pj with < j < p and q^ with 
I < k < q are uniquely determined in terms of the f n 
coefficients with < n < n max by expanding the [p, q] 
Pade approximant in a Taylor series around z = and 
solving the set of n max linear equations. 
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TABLE I: Values of the £-loop beta function coefficients bg defined 
in Eq. J2.2B in the SXJ(N) gauge theory with Ne fermions trans- 
forming according to the fundamental representation, as functions 
of N and Nf , for the range 1 13, II) where the theory is asymptotically 
free. 
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3 


10 


0.345 


-0.156 


-0.386 


0.072 


3 


11 


0.292 


-0.236 


-0.463 


0.154 


3 


12 


0.239 


-0.317 


-0.534 


0.273 


3 


13 


0.186 


-0.397 


-0.599 


0.429 


3 


14 


0.133 


-0.477 


-0.658 


0.622 


3 


15 


0.080 


-0.557 


-0.711 


0.852 


3 


16 


0.0265 


-0.637 


-0.758 


1.121 


4 





1.17 


1.15 


1.71 


3.50 


4 


1 


1.11 


1.04 


1.46 


2.88 


4 


2 


1.06 


0.932 


1.22 


2.31 


4 


3 


1.01 


0.824 


0.986 


1.80 


4 


4 


0.955 


0.716 


0.762 


1.36 


4 


5 


0.902 


0.607 


0.546 


0.972 


4 


6 


0.849 


0.499 


0.339 


0.647 


4 


7 


0.796 


0.391 


0.140 


0.385 


4 


8 


0.743 


0.283 


-0.051 


0.184 


4 


9 


0.690 


0.175 


-0.234 


0.046 


4 


10 


0.637 


0.066 


-0.409 


-0.029 


4 


11 


0.584 


-0.042 


-0.575 


-0.040 


4 


12 


0.531 


-0.150 


-0.733 


0.013 


4 


13 


0.477 


-0.258 


-0.883 


0.131 


4 


14 


0.424 


-0.366 


-1.025 


0.314 


4 


15 


0.371 


-0.474 


-1.16 


0.562 


4 


16 


0.318 


-0.583 


-1.28 


0.877 


4 


17 


0.265 


-0.691 


-1.40 


1.26 


4 


18 


0.212 


-0.799 


-1.51 


1.71 


4 


19 


0.159 


-0.907 


-1.61 


2.22 


4 


20 


0.106 


-1.015 


-1.70 


2.81 


4 


21 


0.053 


-1.124 


-1.79 


3.46 
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TABLE II: Values of Nf M2z , N }fiZZi± , and N fMzJ , i = 2,3, 
for SU(JV) with Nf fermions in the fundamental representation. 
We only list physical, i.e., real, non-negative values. Thus, since 
•N/,6z4,i < 0) i s i s not included. 



N 


^ f,max 


Nf,b2z 


{N f , b 3 Z ,-,Nf, b3z ,+) 


(Nf : b4z,2,Nf : b4z,3) 


2 


ll 


5.55 


(3.99, 27.6) 


none 


3 


16.5 


8.05 


(5.84, 40.6) 


none 


4 


22 


10.61 


(7.73, 53.8) 


(9.51,11.83) 



TABLE III: Estimates of a cr R from the onc-gluon exchange ap- 
proximation to the Dyson-Schwingcr equation for the fermion prop- 
agator. Values are listed for SU(A r ) with 2 < N < 6 and the rep- 
resentations R = (i) fundamental (fund), (ii) adjoint (adj), (iii) 
symmetric rank-2 tensor (S2), and (iv) antisymmetric rank-2 (A2). 



N 


Otcr,fund 


Qtcr.adj 


a cr ,S2 


Ctcr,A2 


2 


1.40 


0.52 


0.52 




3 


0.79 


0.35 


0.31 


0.79 


4 


0.56 


0.26 


0.23 


0.42 


5 


0.44 


0.21 


0.19 


0.29 


6 


0.36 


0.17 


0.16 


0.22 
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TABLE IV: Values of the (approximate or exact) IR zeros in a 
of the SU(iV) beta function with Nf fermions in the fundamental 
representation, for N = 2,3,4, calculated at n-loop order, and 
denoted as cijr n £. For each N, we only give results for the integral 
Nf values in the range l|3,7|l where the theory is asymptotically 
free and the two-loop beta function has an infrared zero. For the 
four-loop beta function, the cubic equation d 3- 32 1 has three zeros, 
one of which is negative, one of which is aiR^i, and the third of 
which is positive but farther from the origin. We include the latter, 
denoted as a^i u . We also list zeros from the [1,2] and [2,1] Pade 
approximants to the four-loop beta function. 



N 


N f 


OtIR,2l 






Q /B,«,[l,2] 


a IRAl,[2,\] 




2 


7 


2.83 


1.05 


1.21 


2.30 


1.16 


4.12 


2 


8 


1.26 


0.688 


0.760 


0.952 


0.741 


3.11 


2 


9 


0.595 


0.418 


0.444 


0.475 


0.438 


2.395 


2 


10 


0.231 


0.196 


0.200 


0.202 


0.200 


1.97 


3 


10 


2.21 


0.764 


0.815 


1.47 


0.807 


5.62 


3 


11 


1.23 


0.578 


0.626 


0.871 


0.616 


3.29 


Q 
O 


1 


U. t 04 


U.4O0 


n A 70 


U.0D1 


U.4DZ 


z.zyo 


3 


13 


0.468 


0.317 


0.337 


0.367 


0.333 


1.78 


3 


14 


0.278 


0.215 


0.224 


0.231 


0.222 


1.48 


3 


15 


0.143 


0.123 


0.126 


0.127 


0.125 


1.29 


3 


16 


0.0416 


0.0397 


0.0398 


0.0398 


0.0398 


1.15 


4 


13 


1.85 


0.604 


0.628 


1.14 


0.625 


6.94 


4 


14 


1.16 


0.489 


0.521 


0.776 


0.516 


3.49 


4 


15 


0.783 


0.397 


0.428 


0.556 


0.422 


2.30 


4 


16 


0.546 


0.320 


0.345 


0.407 


0.340 


1.73 


4 


17 


0.384 


0.254 


0.271 


0.298 


0.267 


1.40 


4 


18 


0.266 


0.194 


0.205 


0.215 


0.203 


1.19 


4 


19 


0.175 


0.140 


0.145 


0.149 


0.145 


1.05 


4 


20 


0.105 


0.091 


0.092 


0.0930 


0.0921 


0.947 


4 


21 


0.0472 


0.044 


0.044 


0.0444 


0.0443 


0.870 



19 



TABLE V: Values of the £-loop coefficients in the series expan- 
sion 1 12. 51 1 for the anomalous dimension j m , as functions of N and 
Nf, for the range 1 13.111 where the theory is asymptotically free. 



N 




Cl 


C2 


C3 


C4 


2 





0.358 


0.318 


0.310 


0.329 


2 


1 


0.358 


0.302 


0.254 


0.234 


2 


2 


0.358 


0.286 


0.195 


0.143 


2 


3 


0.358 


0.270 


0.134 


0.0577 


2 


4 


0.358 


0.254 


0.0712 


-0.0218 


2 


5 


0.358 


0.239 


0.00656 


-0.0952 


2 


6 


0.358 


0.223 


-0.0601 


-0.162 


2 


7 


0.358 


0.207 


-0.129 


-0.222 


2 


8 


0.358 


0.191 


-0.199 


-0.274 


2 


9 


0.358 


0.175 


-0.272 


-0.319 


2 


10 


0.358 


0.1595 


-0.346 


-0.355 


3 





0.637 


0.853 


1.26 


2.03 


3 


1 


0.637 


0.825 


1.11 


1.64 


3 


2 


0.637 


0.796 


0.957 


1.27 


3 


3 


0.637 


0.768 


0.801 


0.909 


3 


4 


0.637 


0.740 


0.642 


0.561 


3 


5 


0.637 


0.712 


0.479 


0.227 


3 


6 


0.637 


0.684 


0.312 


-0.0926 


3 


7 


0.637 


0.656 


0.142 


-0.396 


3 


8 


0.637 


0.628 


-0.0313 


-0.683 


3 


9 


0.637 


0.599 


-0.208 


-0.953 


3 


10 


0.637 


0.571 


-0.389 


-1.21 


3 


11 


0.637 


0.543 


-0.573 


-1.44 


3 


12 


0.637 


0.515 


-0.760 


-1.65 


3 


13 


0.637 


0.487 


-0.951 


-1.85 


3 


14 


0.637 


0.459 


-1.145 


-2.02 


3 


15 


0.637 


0.431 


-1.34 


-2.18 


3 


16 


0.637 


0.402 


-1.54 


-2.31 


4 





0.895 


1.60 


3.17 


6.86 


4 


1 


0.895 


1.56 


2.89 


5.88 


4 


2 


0.895 


1.52 


2.61 


4.93 


4 


3 


0.895 


1.48 


2.32 


4.00 


4 


4 


0.895 


1.44 


2.03 


3.09 


4 


5 


0.895 


1.40 


1.73 


2.21 


4 


6 


0.895 


1.36 


1.43 


1.36 


4 


7 


0.895 


1.33 


1.12 


0.526 


4 


8 


0.895 


1.29 


0.808 


-0.275 


4 


9 


0.895 


1.25 


0.492 


-1.05 


4 


10 


0.895 


1.21 


0.170 


-1.79 


4 


11 


0.895 


1.17 


-0.157 


-2.50 


4 


12 


0.895 


1.13 


-0.488 


-3.18 


4 


13 


0.895 


1.09 


-0.825 


-3.83 


4 


14 


0.895 


1.05 


-1.17 


-4.45 


4 


15 


0.895 


1.01 


-1.51 


-5.025 


4 


16 


0.895 


0.969 


-1.86 


-5.57 


4 


17 


0.895 


0.930 


-2.22 


-6.08 


4 


18 


0.895 


0.890 


-2.58 


-6.54 


4 


19 


0.895 


0.850 


-2.95 


-6.97 


4 


20 


0.895 


0.811 


-3.32 


-7.37 


4 


21 


0.895 


0.771 


-3.69 


-7.71 



20 



TABLE VI: Values of the anomalous dimension in the SU(AT) 
theory with Nf fermions in the fundamental representation, y m , 
calculated to the n-loop order in perturbation theory and evaluated 
at the IR zero of the beta function calculated to this order, cem n £, 
for i = 2,3,4. We denote these as ■y n e(aiR,ni)- For sufficiently 
small Nf > Nf b 2z m each N case, Oijji ^l ls so large that the 
formal value of f2l( Q lR it) is larger than 2 and hence unphysical; 
we indicate this by placing these values in parentheses. 



N 


JV / 




^f3t\ a IR,3£ ) 


y4£\QlR,4£} 


2 


7 


(2 67) 


0.457 


0325 


2 


8 


0.752 


0.272 


0.204 


2 


9 


0.275 


0.161 


0.157 


2 


10 


0.0910 


0.0738 


0.0748 


q 
o 


i n 

-LU 


f4 1 Q"! 

V / 


D R47 


U. LOU 


q 


]^ 


1 fil 


fl 4^Q 


n 9^0 


q 


1 9 


u. it o 


n 9 


U.i JO 


q 
o 




n 4fl4 


n 99n 


n 91 n 


3 


14 


0.212 


0.146 


0.147 


3 


15 


0.0997 


0.0826 


0.0836 


3 


16 


0.0272 


0.0258 


0.0259 


4 


13 


(5.38) 


0.755 


0.192 


4 


14 


(2.45) 


0.552 


0.259 


4 


15 


1.32 


0.420 


0.281 


4 


16 


0.778 


0.325 


0.269 


4 


17 


0.481 


0.251 


0.234 


4 


18 


0.301 


0.189 


0.187 


4 


19 


0.183 


0.134 


0.136 


4 


20 


0.102 


0.0854 


0.0865 


4 


21 


0.0440 


0.0407 


0.0409 



TABLE VII: Values of the (approximate or exact) IR zeros in a 
of the SU(N) beta function with Nf = 2 fermions in the adjoint 
representation, for TV = 2,3,4, calculated at n-loop order, and 
denoted as cxui^^dj- For the four-loop beta function, the cubic 
equation 113.3211 has three zeros, one of which is negative, one of 
which is ctiR,4t !a dj, and the third of which is positive but farther 
from the origin. We include the latter, denoted as 0i4i : u : adQ- We 
also list zeros from the [1,2] and [2,1] Pade approximants to the 
four-loop beta function. 



N 


OLIR,2t,adj 


aiR,3l,adj 


OLIR,4l,adj 


0(IR,4e,[l,2],adj 


a IR, 41, [2,1], adj 


Ot4l,u,adj 


2 


0.628 


0.459 


0.493 


0.522 


0.484 


2.235 


3 


0.419 


0.306 


0.323 


0.341 


0.320 


1.80 


4 


0.314 


0.2295 


0.241 


0.253 


0.239 


1.47 
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TABLE VIII: Values of the anomalous dimension 7 m in an SU(iV) 
gauge theory with Nf = 2 (Dirac) fermions in the adjoint repre- 
sentation, calculated to the n-loop order in perturbation theory 
and evaluated at the IR zero of the beta function calculated to 
this order, for n = 2,3,4. We denote these as "Yne^dji&l R,nl,adj)- 
We also list the value of ^e, adj evaluated at a equal to the /3DS 
estimate, Eq. 1 19.311 , for a criac ij). 



N 


l2l,adj{0!lRfil,adj) 


y3e,adj(0!lR^e,adj) 


7«,ad;(Q!/fl,«,a<ij) 


72f,adj(Q! cr , a[ Jj) 


2 


0.820 


0.543 


0.571 


0.653 


3 


0.820 


0.543 


0.561 


0.653 


4 


0.820 


0.543 


0.557 


0.653 



TABLE IX: Values of the (approximate or exact) IR zero in a of 
the SXJ(N) beta function with Nf = 2 fermions in the symmetric 
rank-2 (i.e., S2) representation, for N = 3,4, calculated at n-loop 
order, and denoted as ct IR n i g2- 



N 


N f 


OtIR,2l,S2 


aiR,3t,S2 


&IR,4t,S2 


3 


2 


0.842 


0.500 


0.522 


3 


3 


0.085 


0.079 


0.080 


4 


2 


0.967 


0.485 


0.485 


4 


3 


0.152 


0.129 


0.132 



TABLE X: Values of the (approximate or exact) IR zero in a of 
the SU(4) beta function with Nf fermions in the antisymmetric 
rank-2 (i.e., A2) representation, for the range 5 < Nf < 10 where 
the theory is asymptotically free and has an IR zero of the beta 
function, calculated at n-loop order, and denoted as otiR in t t A2- 



N 


N f 


OiIR,2l,A2 


OiIR,Zl,A2 


CtlRM,A2 


4 


6 


2.17 


0.664 


0.723 


4 


7 


0.890 


0.437 


0.482 


4 


8 


0.449 


0.287 


0.312 


4 


9 


0.225 


0.174 


0.183 


4 


10 


0.090 


0.080 


0.082 



TABLE XI: Values of 7 m in an SU(A r ) gauge theory with Nf 
fermions in the symmetric rank-2 tensor representation S2, cal- 
culated to the n-loop order in perturbation theory and evaluated 
at the IR zero of the beta function calculated to this order, for 
n = 2,3,4. We denote these as ■y„i S2( a lR nl,S2)- We also list 
72£, S2 evaluated at a equal to the estimate Eq. I|9.3[l for ct cr ,S2- 



N 


N f 


l2t,S2{(XlR,2t,S2) 


l3l,S2{&IR,3,l,S2) 


7«,S2(«/fl,4^,S2) 


72£,S2(«cr,S2) 


3 


2 


(2.44) 


1.28 


1.38 


0.653 


3 


3 


0.144 


0.133 


0.134 


0.619 


4 


2 


(4.82) 


(2.08) 


(2.27) 


0.659 


4 


3 


0.381 


0.313 


0.319 


0.629 
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TABLE XII: Values of r y m in an SU(JV) gauge theory with Nf 
fermions in the antisymmetric rank-2 tensor representation A2, cal- 
culated to the n-loop order in perturbation theory and evaluated 
at the IR zero of the beta function calculated to this order, for 
N = 4 and n = 2,3,4. We denote these as ■y n i,A2(c*l R,nl,A2)- We 
also list 72£,A2 evaluated at a equal to the estimate Eq. 119, 3H for 
a cr ,A2- 



N 


Nf 


J2£.,A2(CXIR,2£,A2) 


~l3l,A2{a-IR,3l,A2) 


r Y4,l,A2\OllR,4l,A2) 


72<?,A2(a C r,A2) 


4 


6 


(9.78) 


1.38 


0.647 


0.769 


4 


7 


(2.19) 


0.695 


0.505 


0.750 


4 


8 


0.802 


0.402 


0.375 


0.732 


4 


9 


0.331 


0.228 


0.231 


0.713 


4 


10 


0.117 


0.101 


0.102 


0.695 
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We consider a vectorial, asymptotically free gauge theory and analyze the effect of higher-loop 
corrections to the beta function on the evolution of the theory from the ultraviolet to the infrared. 
We study the case in which the theory contains Nf copies of a fermion transforming according to 
the fundamental representation and several higher-dimensional representations of the gauge group. 
We also calculate higher-loop values of the anomalous dimension of the mass, j m of {pip at the 
infrared zero of the beta function. We find that for a given theory, the values of 7™, calculated to 
three- and four-loop order, and evaluated at the infrared zero computed to the same order, tend to 
be somewhat smaller than the value calculated to two-loop order. The results are compared with 
recent lattice simulations. 



I. INTRODUCTION 

In this paper we investigate how higher-loop correc- 
tions to the beta function affect the evolution of a vec- 
torial, asymptotically free gauge theory (in (3 + 1) di- 
mensions, at zero temperature) from the ultraviolet to 
the infrared. We assume that the theory contains a cer- 
tain number, Nf, of massless Dirac fermions ip transform- 
ing according to a representation R of the gauge group. 
We consider cases where R is the fundamental, adjoint, 
and rank-2 symmetric or antisymmetric tensor represen- 
tation. We also study the effect of higher-loop correc- 
tions to the anomalous dimension 7 m of the fermion mass. 
This work yields more complete information on the na- 
ture of the evolution of the theory from the ultraviolet to 
the infrared, in particular, on the determination of the 
infrared zero of the beta function and the scaling behav- 
ior of the tptp operator in the vicinity of this zero. We 
will give a number of results for a general gauge group G 
but will focus on the case G = SU(iV). 

We denote the running gauge coupling of the the- 
ory as g(ii), with a(jj,) — g{n) 2 / (47r), where jj, is the 
Euclidean energy /momentum scale (which will often be 
suppressed in the notation). The property that the 
SU(JV) gauge interaction is asymptotically free means 
that lim^oo a(p) = 0, and, since the beta function 
is negative for small a, it follows that, as the en- 
ergy/momentum scale fx decreases from large values, a 
increases. As ll decreases and the theory evolves into 
the infrared, two different types of behavior may occur, 
depending on the fermion content. In a theory with a suf- 
ficiently small number of fermions in small enough rep- 
resentations R, as ii decreases through a scale A, the 
coupling a exceeds a critical value a_R jC r, depending on 
R, for the formation of bilinear fermion condensates, and 
these condensates are produced. This may or may not 
be associated with an infrared (IR) zero of the two-loop 
beta function at a value a = am; if the two-loop (2£) 
beta function does have an IR zero, am,2£, then this 
type of behavior requires that am^e. > aR yCr (H, Q. As 
ll decreases toward A and a increases toward the IR zero 
of the beta function, the increase of a as a function of 



decreasing ll is reduced. This gives rise to an a that is 
of order unity, but varies slowly as a function of ll. This 
behavior is interestingly different from the behavior of 
the gauge coupling in quantum chromodynamics (QCD). 
As the condensates form, the fermions gain dynamical 
masses of order A, so that in the low-energy effective the- 
ory applicable at scales ll < A, they are integrated out, 
and the further evolution of the theory into the infrared 
is controlled by the Nf — beta function. 

Alternatively, if the theory has a sufficiently large num- 
ber Nf of fermions and/or if these fermions are in a large 
enough representation R (as bounded above by the re- 
quirement of asymptotic freedom), then the IR zero of 
the beta function occurs at a value smaller than au tCr , 
so that as the scale ll decreases from large values, the the- 
ory evolves into the infrared without ever spontaneously 
breaking chiral symmetry. In this latter case, the IR zero 
of the beta function is an exact IR fixed point (IRFP), 
approached from below as ll — > 0. More complicated 
behavior occurs in theories containing fermions in sev- 
eral different representations Q; here we restrict to the 
case of fermions in a single representation. For a given 
asymptotically free theory that features an IR fixed point 
at a = am, the value of this IRFP decreases as a func- 
tion of Nf. There is thus a critical value of Nf, denoted 
Nf ;Cr , depending on R, at which am decreases below 
an,cr- This value serves as the boundary, as a function 
of Nf, between the interval of nonzero 1 < Nf < Nf >cr 
where the theory evolves into the infrared in a manner 
that involves fermion condensate formation and associ- 
ated spontaneous chiral symmetry breaking (SxSB), and 
the interval Nf :Cr < Nf < Nf >max , where the theory 
evolves into the infrared without this condensate forma- 
tion and chiral symmetry breaking, with Nf tTnax denoting 
the maximal value of Nf consistent with the requirement 
of asymptotic freedom. 

The anomalous dimension 7 m contains important in- 
formation about the scaling behavior of the operator tpip 
for which m is the coefficient, as probed on different mo- 
mentum scales. In a theory with an am ~ 0(1), it 
follows that j m may also be 0(1), which can produce sig- 
nificant enhancement of dynamically generated fermion 



2 



masses due to the renormalization-group factor 



rj = exp 



7m (<*(//)) 



(1.1) 



In a phase where no dynamical mass is generated, 7 m 
simply describes the scaling behavior of the xpip operator. 

There are several motivations for the study of higher- 
order corrections to this evolution of an asymptotically 
free theory into the infrared. First, the critical coupling, 
Q!_R,cr, is generically of order unity, and hence there is a 
need to have a quantitative assessment of the importance 
of higher-loop corrections to the evolution of the the- 
ory. Second, besides being of fundamental held-theoretic 
interest, a knowledge of this evolution plays an impor- 
tant role in modern technicolor (TC) models with dy- 
namical electroweak symmetry breaking, in which the 
slow running of the coupling associated with an approx- 
imate infrared zero of the beta function provides nec- 
essary enhancement of quark and lepton masses [l|, [|[ 
(recent reviews include [4|- @), and can reduce techni- 
color corrections to precision electroweak quantities 
In addition to the fundamental representation, fermions 
in higher-dimensional representations have been studied 
in the context of technicolor [f| Q . Fermions in higher- 
dimensional representations of chiral gauge groups have 
long played a valuable role in studies of extended tech- 
nicolor (ETC) models that were reasonably ultraviolet- 
complete and explicitly worked out the details of the se- 
quential breaking of the ETC chiral gauge symmetries 
down to the TC group 10]. Recently, there has been 
a considerable amount of effort devoted to lattice stud- 
ies of gauge coupling evolution and condensate forma- 
tion in vectorial SU(iV) gauge theories as a function of 
Nj^Jot fermions in both the fundamental representation 
(8|,[l]J-[lj|and higher representations [ll|-[23[ (a recent 
review is [24]). Thus, another important motivation for 
the present work is to provide higher-order calculations 
that can be compared with these lattice studies. 



II. GENERAL THEORETICAL FRAMEWORK 



A. Beta Function 



series for other quantities in quantum field theories do not 
have finite radii of convergence but are only asymptotic, 
experience shows that in situations where the effective ex- 
pansion parameter (here, (ot/n) times appropriate group 
invariants) is not too large, the first few terms can pro- 
vide both qualitative and quantitative insight into the 
physics. The first two coefficients in the expansion (12.21) . 
which are scheme- independent, are [25j 

h = hllC A -4T f N f ) (2.3) 



and |'2f 



&2 = - [3ACl-A(5C A + 3C f )T f N f ] . (2.4) 

Here C/ = C2(R) is the quadratic Casimir invariant for 
the representation R to which the Nf fermions belong, 
Ca = C*2(G) is the quadratic Casimiar invariant for the 
adjoint representation, and Tf = T(R) is the trace in- 
variant for the fermion representation R. Higher-order 
coefficients, which are scheme-dependent [13]) have been 
calculated up to four-loop order [H, |29| _. Some further 
details are given in Appendix I. Values of bi fof 1 < TV < 4 
and relevant rangles of Nf are given in Table [H 



B. Anomalous Dimension of the tpip Operator 

The anomalous dimension 7 m for the fermion bilinear 
iptp describes the scaling properties of this operator and 
can be expressed as a series in a or equivalently, a, 



c e a 



(2.5) 



where ct — ci/(An) 1 is the ^-loop series coefficient. Via 
Eq. (|1.1[) , the anomalous dimension 7 m governs the run- 
ning of a dynamically generated fermion mass. The co- 
efficients C£ have been calculated to four-loop order [30l ]. 
The first two are 



The beta function of the theory is denoted f3 = dg/dt, 
where dt = dln/i. In terms of the variable 



9 



16tt 2 4tt 



(2.1) 



the beta function can be written equivalently as /3 a = 
da/dt, expressed as a series 



d 

— = -2a be a e — -2a ^^b^a 1 , 



(2.2) 



where I denotes the number of loops involved in the cal- 
culation of bi and b( — bi/(A-K) 1 . Although this series and 



and 



c 2 = 2C f 



ci = 6C f 



3 r + 97 r 10 t n 
-C f + -C A - -T f N f 



(2.6) 



(2.7) 



For reference, the coefficient C3 is listed in Appendix I. 
Since as V/ approaches Nf, max from below, b\ — > with 
nonzero 62 and hence aiR — > 0, and since the perturba- 
tive calculation expresses 7 m in a power series in a, it 
follows that as 7 m — > as Nf approaches Nf_ max from 
below. We note that a conjectured beta function that 
directly relates to 7 has been proposed (3l| . 
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III. PROPERTIES OF BETA FUNCTION 
COEFFICIENTS AND APPLICATION TO 
FUNDAMENTAL REPRESENTATION 

In this section we discuss some general properties of 
the beta function coefficients as functions of Nf, and 
give particular results for the case of fermions in the fun- 
damental representation. In later sections we consider 
fermions in two-index representations. 



61 



Since we restrict our considerations to an asymptoti- 
cally free theory, we require that, with our sign conven- 
tions, b\ > 0. This, in turn, implies that 



N f < N f , 



where 



Nf, 



UCa 
4T f 



(3.1) 



(3.2) 



Thus, for fermions in the fundamental representation, 

Nf,max,fund = (11/2)JV. 



B. 62 and Condition for Infrared Zero of (3 

We next proceed to characterize the behavior of the 
higher-loop coefficients of the beta function, bi with 
I = 2,3,4, and the resultant zero(s) of the beta func- 
tion, in terms of their dependence on Nf. The two- loop 
results are well-known and are included here so that the 
discussion will be self-contained. Since only the first two 
coefficients of the beta function are scheme- independent, 
it follows that, to the extent that one is in a momentum 
regime where one can reliably use the perturbative beta 
function, the zeros obtained from these first two coeffi- 
cients should be sufficient to characterize the physics at 
least qualitatively. When one includes higher-loop con- 
tributions to the beta function, one expects shifts of ze- 
ros, and there are, indeed, generically substantial shifts 
if zeros of the two-loop beta function occur at a ~ O(l). 
However, if inclusion of three- and/or higher-loop contri- 
butions to P leads to a qualitative change in behavior, 
relative to the behavior obtained from the two-loop /3 
function, then the results cannot be considered fully reli- 
able, since they are scheme-dependent. For example, for 
a given gauge group G and fermion content, if the two- 
loop beta function does not have an infrared zero but 
the three-loop beta function does, one could not con- 
clude reliably that this is a physical prediction of the 
theory. Moreover, it should be noted that even if there is 
no zero of the two-loop beta function away from the ori- 
gin, i.e., a perturbative IRFP, the beta function may ex- 
hibit a nonperturbative slowing of the running associated 
with the fact that at energy scales below the confinement 



scale, the physics is not accurately described in terms of 
the Lagrangian degrees of freedom (fermions and gluons) 

[al-HS. 

Another general comment is that the expression of the 
beta function in Eq. (|2 . 2[) is semiperturbative and does 
not incorporate certain nonperturbative properties of the 
physics, such as instantons, whose contributions involve 
essential zeros of the form exp(— kit /a), where k is a nu- 
merical constant. These instanton effects are absent to 
any order of the perturbative expansion in Eq. (|2.2p but 
play an important role in the theory. For example, they 
break the global U(1)a symmetry [35| and also enhance 
spontaneous chiral symmetry breaking (36[-[3!|. Esti- 
mates of the effects of instantons on the running of a in 
quantum chromodynamics have found that they increase 
this running, i.e., they make j3 more negative in the re- 
gion of small to moderate a values [37|- If one were to 
model the effect of instantons crudely via a modification 
of j3 such as 



_da 
Pa ~ dt 



-2a' 



^6,a £ - 1 +Aexp(-^) 



(3.3) 



then, since A > 0, this would have the effect of increasing 
the value of the smallest (nonzero, positive) IR zero am 
of p. For a given minimal value of a cr .R for condensate 
formation and spontaneous chiral symmetry breaking, 
since at least at the perturbative level am is a decreas- 
ing function of Nf, it would follow that incorporating 
instanton effects would increase the value of Nf iCr , i.e., 
would increase the interval in Nf where there is SxSB. 
Furthermore, since instantons enhance chiral symmetry 
breaking, they would tend to reduce the value of a cr .R, 
which also has the same effect of increasing Nf tCr . We 
shall comment below on how, although the semipertur- 
bative one-gluon exchange approximation to the Dyson- 
Schwinger (DS) equation does not directly include effects 
of confinement or instantons, it may nevertheless yield an 
approximately correct value of Nf jCr because of another 
approximation involved that has the opposite effect on 
the estimate. 

If one knows the beta function calculated to a maximal 
loop order £ max , then the equation for the zeros of the 
beta function, aside from the zero at a = 0, is 



E 6 

i=i 



= b l 



1 + 



1=1 



. (3.4) 



As is clear from Eq. (|3.4[) . the zeros of j3 away from 
the origin depend only on the £ ma x — 1 ratios bifb\ for 

^ — ^ — @max ■ 

The coefficients b\ and 62 are linear functions of Nf, 
while 63 and 64 are, respectively, quadratic and cubic 
functions of Nf. With our sign convention in which an 
overall minus sign is extracted in Eq. (|2.2p , each of these 
coefficients is positive for Nf = 0. The coefficients &i and 
62 are both monotonically and linearly decreasing func- 
tions of Nf. As Nf increases sufficiently, 62 thus reverses 



4 



sign, from positive to negative, vanishing at Nf = Nf^2zi 
where 



N 



f,b1z 



2T f (5C A + 3C f ) 



(3.5) 



(The subscript biz stands for the condition that bg is 
zero). Since 



N 



f,max 



N 



}.b2z 



ZC A {\1C S + 7C A ) 
42>(3C> + 5CU) 



> 



(3.6) 



i.e., Nf^ max > Nf,b2z, it follows that there is always a 
nonvacuous interval in the variable Nf where the theory 
is asymptotically free and the two-loop (2i) beta function 
has an infrared zero, namely 



N f , b2z <N f < N f , 



This zero occurs at 



OLIR.21 



A-nbi 



(3.7) 



(3.8) 



and is physical for b 2 < 0. Explicitly, for the fundamental 
representation, 



Nffi2zjund 



UN A 



UN 2 - 3 



and 



&IR,2e,fund 



4ir{llN-2N f ) 
-347V 2 + N f (lSN - 37V- 1 ) 



(3.9) 



(3.10) 



Illustrative values of Nf^2z,fund are given in Table ITT] 
The sizes of ^-loop contributions arc determined by 
(a /it) 1 multiplied by corresponding powers of various 
group invariants. Illustrative values of uiR^njund are 
given in Table [TTT] for N = 2, 3, 4 and the subset of the 
interval (13.71) for which am^ejund is not so large as to 
render the two-loop perturbative calculation obviously 
unreliable. Here and below, when a and 7 values are 
listed without an explicit R, it is understood that they re- 
fer to the fundamental representation. Examples of cases 
that we do not include in the table because the two-loop 
result cannot be considered reliable include the follow- 
ing (with formal values of am^ejund listed): N = 2, 
N f = 5, where ai R:2 ejund = H-4; N = 3, Nf = 9, 
where am^ejund — 5.2; and N = 4, Nf = 11, 12, where 
0!IR,2£Jund — 14, 3.5. 



For reference, the estimate in Eq. (J9T3J) of a cr from the 
analysis of the Dyson-Schwinger equation for the fermion 
propagator, in the one-gluon exchange approximation, 
takes the form in Eq. (|9.3[) for a fermion in the funda- 
mental representation. As listed in Table HVl this has the 
respective values 1.4, 0.79, and 0.56 for N = 2, 3, 4, re- 
spectively (where we quote the results to two significant 
figures but do not mean to imply that they have such 
a high degree of accuracy). Setting a cr = aiR,2t yields 
the resultant estimates of Nf^ cr , which, rounded to the 
nearest integers, are 8, 12, and 16 for these values of N. 
We denote these as (3DS estimates since they combine 
a calculation of ajR from the perturbative two-loop j5 
function with the (one-gluon exchange approximation to 
the) Dyson-Schwinger equation. 



As Nf approaches its maximum value, Nf tmax , al- 
lowed by the constraint that the theory be asymptot- 
ically free, 62 reaches its most negative value, namely 
b 2 = —C A (7C A + HC/). Clearly, for Nf values such 
that 62 is only negative by a small amount and am t 2t 
is large, the perturbative calculation is not reliable. As 
Nf increases further in the range (|3.ip and uir^i de- 
creases, the calculation becomes more reliable. In Table 
HI1 we list the numerical values of Nf^2z for some illus- 
trative values of N . At the two- loop level, depending on 
whether azR,2i is smaller or larger than a critical value 
for fermion condensation, this is an exact or approxi- 
mate infrared fixed point (IRFP) of the renormalization 
group for the gauge coupling. The existence of such an 
IRFP is of fundamental importance in determining how 
the theory evolves from the ultraviolet to the infrared 
[40I ]. In particular, as mentioned above, this determines 
whether, as the scale \i decreases sufficiently to a scale 
A (depending on the group G and the fermion content), 
a grows to a large enough size to produce fermion con- 
densates or, on the contrary, the coupling never gets this 
large and the theory evolves into the infrared in a chirally 
symmetric manner, without ever producing such fermion 
condensates. Note that in the former case, the fermions 
involved in the condensates get dynamical masses of or- 
der A and are integrated out of the effective low-energy 
field theory applicable for scales /i < A, so that the fur- 
ther evolution into the infrared is governed by a different 
beta function. 



It is useful to observe how rapidly the numbers Nf^2z 
approach their large-A values. The number Nf^zjund 
has the large- iV expansion 



Nf : b2z,fund = N 



34 
13 



102 



306 



(13A) 2 (13) 3 A^ 4 



O 



N 6 J 



= N 



2.615 



0.60355 0.1393 



N 2 



N 4 



O 



N e J 



(3.11) 



As is evident from Table [TH the values of Nf^2zjund approach the leading asymptotic form for moderate A, as a 
result of the fact that the subleading term in Eq. (|3.11[) is suppressed by 1/N 2 . 
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It is of interest to consider the 't Hooft large- TV limit, where 

N -> oo with aN fixed. (3.12) 

In a theory with fcrmions in the fundamental representation, in order for them to have a non-negligible effect in this 
limit, one considers the simultaneous Veneziano limit 

N f -> oo with r=^- fixed. (3.13) 

In the combined limit of Eqs. (|3.12|) and (|3. 131) . the range of r satisfying the requirement of asymptotic freedom and 
the condition that hi < so that the two-loop beta function has an IR zero is [41j 

^ < r < — ' , i.e., 2.615 <r< 5.5. (3.14) 



C. Coefficient bs and Three-Loop Behavior of the Beta Function 

The three-loop beta function coefficient 63 is a quadratic function of Nf with positive coefficients of its iV2 and N'j 
terms and a negative coefficient of its Nf term. Hence, regarded as a function of the formal real variable Nf, it is 
positive for large negative and positive Nf, and positive at Nf = 0. The derivative of b 3 with respect to Nf is 



db 3 
dN f 



1415 
~27~ 



C 



205 



C A C { 



■2Cj 



T 



38 316 
jC f + — C A 



(3.15) 



For the fermion representations R that we consider here, for small values of Nf, this derivative db^/dNf is negative, 
so that in this region of Nf, 63 decreases from its positive value at Nf = as Nf increases. Because 63 is a quadratic 
polynomial in Nf, the condition that it vanishes gives two formal solutions for Nf, namely 



N 



{U15C 2 A 



f,b3z,± 



615C A C f - 5AC] ± 377^3" ) 



47/(7904 + 66C/) 



(3.16) 



and 



F Rb3 = \Tl\hlC\ + 109578C^C*/ + 25045CiC : 



7380CUC? 



3240 



(3.17) 



Given that Fnb3 > 0, as is the case here, so that the 
values Nf t b3 Z j are real, it follows that b 3 is positive in 
the intervals Nf < Nf^ b3z ^_ and Nf > Nf^ b3z ^ + and neg- 
ative in the interval Nf^z,- < Nf < Nf ba z +. The 
value Nf ib 3 Zi+ and the neighborhood of Nf values in the 
vicinity of Nf tb3Zt+ are not of interest here because they 
are larger than the maximal value Nf^ max allowed by the 
requirement of asymptotic freedom, 



Nf,b3z,+ > Nf,max ■ 



(3.18) 



Thus, 63 only changes sign once for Nf in the asymptot- 
ically free interval < Nf < Nf max . As Nf approaches 
Nf m ax from below, 63 decreases to a negative value given 
by' 



\bs)N f =Nf. 



Ca 
24 



1127C1 



UCf(UC A -SC f ) 

(3-19) 

For fermions in the fundamental representation, this is 



{ba)N f =Nf,- 



701 
12 



N" 



121 
12 



-N r . 



11 

8N~ C 



(3.20) 



As is clear from Table UH for this case 

Nf,b3z,l < ^f,b2z ■ 



(3.21) 



We noted above that any physically reliable zero of the 
beta function must be present already at the level of the 
two-loop beta function, since this is the maximal scheme- 
independent part of this function. Hence, in analyz- 
ing such a zero for the case under consideration where 
the fermions transform according to the fundamental 
representation of SU(iV), we only consider the interval 
(|3.7[) . Combining this fact with our results (13.211) and 
(|3.18p . it follows that 63 is negative throughout all of 
the interval (|3.7j) of interest here. For this fundamental- 
representation case, the Nf^3 Z ,j with j = 1, 2 have the 
large- N expansions 



N 



f,b3z,l 



= N 



1.911 + 



0.3244 0.06844 



N 2 



N 4 



O 



ye) 



(3.22) 
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and 

Nf,b3z,2 



N 



13.348- 



1.667 0.3978 



N 2 



N 4 



O 



ye) 



(3.23) 



Here, N f , max = 5.5N. 

At three-loop order, the equation for the zeros of the 
beta function, aside from the zero at a = 0, is b\ + b 2 a + 
63 a 2 = 0. Formally, this equation has two solutions for a 
and hence for a, namely 



2tt 

at/3z,3e,± = -7— 
03 



-bn± 



4&l63 



(3.24) 



Since 62 must be negative in order for the beta function 
to have a scheme-independent infrared zero, and since 
for fermions in the fundamental representation we have 
shown that 63 < in the relevant interval Q3.7p . we can 
rewrite this equivalently as 



2tt 

N 



(3.25) 



In order for a given solution to be physical, it must 
be real and positive. As is evident from Eq. (|3.25j) . the 
solution corresponding to the the + sign in Eq. (|3.24j) 
(i.e., the — sign in Eq. p.25[) ) is negative and hence 
unphysical. Thus, there is a unique physical solution 
for the IR zero of the beta function to three-loop order, 
namely 



&IR,3l = o^ z ,3£, 



(3.26) 



Illustrative values for this IR zero of the beta function at 
three-loop order are listed in Table IIIII 



For an arbitrary fermion representation for which (3 
has a two- loop IR zero, we observe that the value of this 
zero decreases when one calculates it to three-loop order, 
i.e., 



This can be proved as follows. We have 
2tt 



(3.27) 



IM3 



2h\b 3 



bl-MxIbl + MAbs 



(3.28) 

The expression in square brackets is positive if and only 
if 



(2b 1 \b 3 \+b 2 ) 2 -b 2 2 (^+^b 1 \b 3 \) >0 



(3.29) 



But the difference in Q3.29[) is equal to the positive- 
definite quantity 6f 6§ , which proves the inequality (|3.27p . 
This inequality is evident in Tabic IIIII 



D. Coefficient 64 and Four-Loop Behavior of (5 

The four-loop beta-function coefficient, 64, was calcu- 
lated in Ref. [20|. We next analyze its behavior as a 
function of Nf and the result four-loop IR zero of the 
beta function. The coefficient 64 is a cubic polynomial 
in Nf which has positive coefficients of its iV°, and N^ 
terms. Hence, regarded as a function of the formal real 
variable Nf, 64 is negative for large negative Nf, positive 
for Nf — 0, and also positive for large positive Nf. For 
fermions in the fundamental representation, the deriva- 
tive at Nf = is 



P-) = (™™ + ^(3)) N* + ( ^ - ^ C (3)) N + (- 2 -^ + ^C®) AT" 1 - ^A- 3 (3-30) 
dN f J n _ o \ 1944 9 ^ v 7 V !944 9 SV ') \ 216 9 ^ 7 8 v ' 



where C(z) is the Riemann zeta function, 

00 

n— 1 

and C(3) = 1.20205690.... This derivative is negative for 
all N . (In the complex N plane, it has six zeros at three 
complex-conjugate pairs of N values.) It follows that, 
again as a function of the formal real variable Nf, 64 
has a local maximum at a negative value of Nf and then 
decreases through positive values as Nf increases toward 
and passes through into the interval of physical values. 

The detailed behavior of &4 in the physical asymp- 
totically free interval < Nt < Nf :Tnax depends on 



N . In particular, one may determine the value of Nf 
where 64 has a minimum and whether 64 has any ze- 
ros for positive Nf. For SU(2), 64 decreases to a mini- 
mum positive value as Nt ascends through the approxi- 
mate value Nf = 5.8, and then increases monotonically 
for larger Nf, so that it is positive-definite for all non- 
negative Nf, in particular, the asymptotically free region 
< Nf < 11. For SU(3), 64 is also positive-definite 
for all non- negative Nf, reaching a local minimum as 
Nf ascends through a value of approximately 8.2 and 
then increasing monotonically for larger Nf. However, 
for SU(4), 64 is positive for < Nt < 9.51, negative for 
the interval 9.51 < Nf < 11.83, and positive again for 
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Nf > 11.83, with zeros at N f ~ 9.51 and N f ~ 11.83. 
We list these zeros of 6 4 as a function of iVf in Table [TTJ 
(Again, we recall that the physical values of Nf are, of 
course, restricted to non-negative integers.) Thus, this 
reversal of sign occurs in the interval of interest here, 
< Nf < 22, where the SU(4) theory is asymptotically 
free. For SU(5), 64 behaves in a manner qualitatively sim- 
ilar to the SU(4) case; it is positive for < Nf < 11.18, 
negative in the interval 11.18 < Nf < 15.18, and posi- 
tive for larger values of Nf, vanishing at Nf ~ 11.18 and 
Nf ~ 15.18. Thus, again, 64 reverses sign in the region 
< Nf < 22.5 where the SU(5) theory is asymptotically 
free. Thus, in contrast with 62 and 63, which are negative 
throughout the interval of Nf of interest (and 61 , which 
is positive), 64 can, for N > 4, vanish and reverse sign in 
this interval. 

At the four-loop level, the equation for the zeros of the 
beta function, aside from a = 0, is the cubic equation 



bi + b 2 a + 6 3 a 2 + b 4 a 3 = 



(3.32) 



This equation has three solutions for a and hence for 
a, which will be denoted ap Z Ai,j, j — 1,2,3. Since the 
coefficients be are real, there are two generic possibilities 
for these three roots, namely that they are all real, or that 
one is real and the other two form a complex-conjugate 
pair. The properties of the roots are further restricted 
by the asymptotic freedom condition that b± > 0, the 
existence of a two-loop IR zero, which requires that 62 < 
0, and the fact that, as we have shown, for the relevant 
range (|3.7p of Nf, where these conditions are met, 63 < 0. 
As is evident in Table IIII1 we find that for the values of 
N and Nf that we consider, the roots of Eq. (|3.32p arc 
real. For all of the values of N and Nf where there is a 
reliable two-loop value for an IR zero of the beta function 
(i.e., where it does not occur at such a large value of a 
as to render the perturbative calculation untrustworthy) , 
one of these roots is negative and hence not physical, one 
of them, namely the minimal positive one, is the physical 
IR zero, which we will denote a,m,4Jl — am 44 /(An), and 
there is a third root at a larger positive value. This third 
root, denoted a^e.u = «4£,«/(47r), is not relevant for our 
analysis, since the initial value of a at a high energy scale 
\x is assumed to be close to zero, so that as the scale fj, 
decreases, a increases and approaches the (positive) zero 
of the beta function closest to the origin, namely am u 
421. 



It is straightforward to display the analytic expressions 
for the root am^i, but we shall not need this for our 
analysis. We list numerical values for am t n for various 
values of N and Nf in Table [TTT1 For completeness, we 
note the specific sets (N,Nf) where am,ii is so large 
that we consider the analysis via the perturbative beta 
function unreliable: these are (N,Nf) = (2,6), (3,9), 
(4,11), and (4,12). 

E. Estimates of Zeros of the Four-Loop Beta 
Function via Pade Approximants 



For the beta function, or more conveniently, 
the reduced function with the prefactor removed, 
Sj™^ -1 bja? -1 , it is useful to calculate and analyze 
Pade approximants, since these provide closed-form ex- 
pressions that, by construction, agree with the series to 
the maximal order to which it is calculated. The ex- 
pansion for (3 a to £ = 4 loop order can be used in two 
ways. First, one can simply solve the cubic equation 
j3 a = bi + b 2 a+b 3 ,a 2 +b^a i = and obtain the three roots, 
one of which is the root of interest, giving the IR zero. 
Secondly, one can calculate Pade approximants, e.g., the 
[2,1] and [1,2] approximants, and determine their zeros. 
The [1,2] Pade approximant has a single zero at 



a Pz,M, [1,2] 



am t u,[i,2] _ 61(6463 - b\ 



4?r b\ - 2616263 + b\b A 



(3.33) 



Taking into account the fact that 62 and 63 are negative 
in the relevant interval (I3.7[) . this can be rewritten as 



a Pz,4l,[l,2\ 



a /ii,«,[l,2] 



MbxM+bl) 



At: |6 2 | 3 + 26 1 |62||6 3 |-6264) 



The two zeros from the [2,1] approximant are 



(3.34) 



,U,[2,1],± 



6263 - 6164 ± 



(6263 -6i6 4 ) 2 - 46x63(62-6264) 



1/2 



2(6264 - 62) 



(3.35) 



Taking account of the fact that 62 and 63 are negative in 
the relevant interval (|3.7|) . this can be rewritten as 



a /3z,«,[2,l],± 



M 4 - NN T (NN - 6iM 2 + 4641631(61 + |6 2 |6 4 ) 



1/2 



2(|6 2 |6 4 
I 



(3.36) 



The expression in Eq. (|3.36|) with the — sign in front expression with the -I- sign in front of the square root 
of the square root is negative and unphysical, while the yields the estimate of the IR fixed point, as ajji^n,!] — 



8 



47ra j a aj 4<,[2,il' As is evident from Eqs. (|3.33[) and (|3.35l) . 
the zeros of the [1,2] and [2,1] Pade approximants incor- 
porate information on /3 up to four loops. One readily 
verifies that in the limit 64 — > 0, the zero of the [1,2] Pade 
reduces to the two- loop result a = —61/62, and the two 
zeros of the [2, 1] Pade reduce to those obtained from the 
three- loop beta function, (|3.26[) . We list the values of 
am obtained from the zeros of the [1,2] and [2, 1] Pade 
approximants to the four-loop beta function for the case 
of fermions in the fundamental representation in Table 

eh 

From our calculations of am at the three- and four- 
loop level for SU(AT) with fermions in the fundamen- 
tal representation, we can make several remarks. Al- 
though n-loop calculations of the beta function for 11 > 3 
loops are scheme-dependent, the results obtained with 
the present MS scheme provide a quantitative measure 
of the accuracy of the scheme-independent two-loop re- 
sult. For a given N, as Nf increases above the mini- 
mal value Nf^2z, where the IR zero first appears, and as 
the resultant ctmu decreases to values ^ 1, the differ- 
ence between <xm,ii and the higher-loop values am,nt f° r 
n = 3, 4 decrease. As is evident from Table Hill the value 
of 

ctiR,n£ generically decreases as one goes from n — 2 to 
n = 3 loops and then increases by a smaller amount as 
one goes from n — 3 to n = 4 loops, so that amM is 
smaller than am,2i- In the same region of Nf values 
such that am t 2£ is reasonably small, the values obtained 
via the [1,2] and [2,1] Pade approximants to the four-loop 
beta function are close to those obtained from the zeros 
of this beta function itself. 



IV. EVALUATION OF THE ANOMALOUS 
DIMENSION 7m AT THE INFRARED ZERO OF /3 

In this section we evaluate the anomalous dimension of 
7 = 7 m , calculated to the n-loop order in perturbation 

I 



We list numerical values of 7(07/?, 2^) in Table IVTl for the 
illustrative values N — 2, 3, 4 and, for each N, a set of 
Nf values in the range Q3.7p . For sufficiently small Nf > 
Nf,b2z in each N case, am,U is so large that the formal 
value of j2i(cyiR,2e) is larger than 2 and hence unphysical; 
we enclose these values in parentheses to indicate that 



theory, at the (approximate or exact) IR zero of the beta 
function to this order, am, n i, for n = 2, 3, 4. We denote 
these as j n e(oiiR t ne)- We focus here on general results 
and their application to the case of fermions in the fun- 
damental representation, and discuss higher-dimensional 
representations in subsequent sections. In general, this 
anomalous dimension must be positive to avoid unphys- 
ical singularities in fermion correlation functions. The 
coefficients ci that enter in Eq. (|2.5I) used in this calcu- 
lation are listed in Table [V] 

A running fermion mass, S(fc), that is dynamically 
generated at a scale A, decays with Euclidean momen- 
tum k > A like 

£(fc)~A(-J (4.1) 

up to logs. Since for k > A, the running coupling a 
is smaller than the critical value a_R jCr and there is no 
spontaneous chiral symmetry breaking, it follows that 
must decrease toward zero as fc/A — > 00. In turn, 
this implies that j m < 2. Hence, a physical value of 7 m 
must lie in the range 

< 7 m < 2 . (4.2) 

For values of Nf such that the theory evolves into the 
infrared in a chirally symmetric manner, so that the IR 
zero of the beta function is exact, the same upper bound 
follows from a related unitarity consideration [431 ] . 

Using the two-loop result for 7 and evaluating it at 
the two-loop value of the IR zero of the beta function, 
we have 



(4.3) 



(4.4) 



they are unphysical artifacts of a perturbative calculation 
at an exessively large value of a. 



In the large- N, large- N f limit of Eqs. (j3"T2l and (|3~T5)) 



C S {\\C A - 4T / A / )(455Cl + 99C A C f + (180(7/ - 2i8C A )T f N f + 80TfNf) 



\2(-nC\ + {\QC A + QCf)T f NfY 
For the fundamental representation, this is 

(N 2 - 1)(1LV - 2N f )(1009N 3 - 99N - (158N 2 + 90)N f + AONNj) 



l2i{aiR,2i) 



12(-34iV 3 + (13iV 2 - 3)AT/) 2 
I 
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with r = Nf/N, Eq. (|4.4I) reduces to 



(11 - 2r)(1009 - 158r + 40r 2 ) 
12(-34+ 13r) 2 



+o 



1 

TV 2 
(4.5) 

For r = 4 corresponding to the asymptotic value of 
N f , crJmd in Eq. (|93I), TSM(ajfl,a<) = 113/144 ~ 0.785, 
which is the same as the large- AT limit of Eq. (|4.7|) . 

One may evaluate rs j2i{p-iR,ii) at ATy equal to the value 
Nf^crjund predicted by the one-gluon exchange (ladder) 
approximation to the Dyson-Schwinger equation for the 
fermion propagator, given in Eq. (|9.4|) . This is somewhat 
formal, since these values of Nf iCr j un d are not, in general, 
integers and hence not actually physical; for example, 
N Ur j md = 7.86, 11.91, 15.94 for N = 2, 3, 4). This 
procedure yields the result 



l2l(o>m^2f, Nf,cr,furvi) 



565iV 4 - 7067V 2 + 225 
1U(N 2 - l)(5iV 2 - 3) 



(4.6) 



For the illustrative cases N — 2, 3, 4, this anomalous 
dimension takes the values 0.88, 0.82, and 0.80, respec- 
tively. As N — > oo, Eq. (|4.6[) has the expansion 



113 

72Ha C r,/ t md) = Y44 



11 



40 A^ 2 



(4.7) 



Since the estimate (|9.4p is close to 4AT even for the small- 
est value, N = 2, and asymptotically approaches AN as 
A" — > 00, it is worthwhile to compare the above values 
of 7, viz., 0.88, 0.82, and 0.80 for AT = 2, 3, 4, with 
72£(a/fl,2£) evaluated at the nearest physical, integer val- 
ues of Nf, namely N f = 8, 12, 16 for N = 2, 3, 4. 
This procedure yields ^iit{&iR,2i) = 0.75, 0.77, 0.78, as 
recorded in Table [VTl To within the strong-coupling the- 
oretical uncertainties of these calculations, these values 
are mutually consistent. 

A closely related approach is to evaluate the two-loop 
expression for j m at a = a crt R, where a cr .n is the 
estimate of the critical coupling for fermion condensa- 
tion obtained from the one-gluon exchange approxima- 
tion to the Dyson-Schwinger equation, and then substi- 
tute Nf — Nf_ cr from the /3DS analysis (see Appendix 
II). This yields the result 



lii{a cri R; N f = Nf jCr , R ) 



21C\ + 128C A C'f + 225C 



J 



lUCf(C A + 3C f 



(4.8) 

For the fundamental representation, this reduces to the 
same result as was obtained in Eq. (|4.6p . 

We have evaluated the three-loop result for 7 at the 
three-loop value of the IR zero of the beta function, which 
we denote as 73^ {olir, 31), and the four-loop result for 7 
at the four-loop value of the IRFP, which we denote as 
r yu{<xiR,u)- We list the resultant values in Table I VII 
From our calculations of 7 m for the case of fermions 
in the fundamental representation, we can make sev- 
eral observations. Although computations of aiR^ni and 
Jne{ a IR,ni) are scheme-dependent for n > 3 loops, they 



provide a useful measure of the accuracy of the lowest- 
order results. As was the case with the position of Q.iB,,nt 
itself, we find that, for a given N and for Nf reasonably 
well above Nf^z so that the perturbative calculation of 
a iR,nt is not too large, the value of Jne( a iR,ni) generi- 
cally decreases as one goes from n = 2 to n = 3 loops. 
Some of this decrease can be ascribed to the decrease in 
&iR,ni going from (n = 2)-loop to (n = 3)-loop order. At 
the four-loop level, ju{atiR,4t) tends to be smaller than 
l3t(®iR,3t) f° r values of Nf from Nf^ 2z to values of Nf 
slightly above the middle of the range Q3.7p . while for 
values of Nf in the upper end of this range, jui&iR^i) 
is slightly larger than 73^(077?, 3^). In general, for the 
values of Nf where am is sufficiently small that the cal- 
culation may be trustworthy, the value of the anomalous 
dimension evaluated at the IR zero of the beta function 
(both calculated to n-loop order) ynt{ a IR,nt)i is some- 
what smaller than unity. 

Several recent high-statistics lattice simulations have 
been carried out on an SU(3) gauge theory with a vary- 
ing number Nf of fermions in the fundamental represen- 
tation in the range 6 < N f < 12 [Il|-[ll, [|, ||. This 
work has yielded evidence for a regime of slowly running 
gauge couplings for Nf ^ 12, consistent with the pres- 
ence of an IR zero of the beta function, in agreement 
with the earlier continuum estimates in Ref. 1]. Ref. 
[ill ] also found a considerable enhancement of (i/'V')//p 
in the SU(3) theory with Nf — 6. Further lattice sim- 
ulations and analysis of data should yield values of j m 
that can be compared with our higher-loop calculations 
in this paper. A preliminary study of the SU(2) theory 
with Nf = 6 fermions has also been reported [2l|. 

V. ADJOINT REPRESENTATION 

In this section we analyze the SU(A r ) theory with Nf 
copies of a Dirac fermion, or equivalently, 2Nf copies of 
a Majorana fermion, in the adjoint representation. For 
this case, the general expression for the maximal value 
of A^ allowed by the requirement of asymptotic freedom, 
Eq. (f3~2")) , reduces to 



N 



f \max ,ad 



11 
T 



(5.1) 



i.e., restricting Nf to the integers, Nf^ max = 2. The gen- 
eral expression in Eq. (|3.5[) for the value of Nf at which 
62 changes sign from positive to negative with increasing 
Nf reduces to 



N ftb2Ztadj = ^ = 1.0625 



(5.2) 



Hence there is only one (integer) value of Nf, namely 
2 Dirac fermions (equivalently, Nf = 4 Majorana 



Nf 

fermions), for which the theory is asymptotically free and 
has an IR zero of the two-loop beta function. This zero 
occurs at 

2-k 1 257 

0tIR,2l,adj = 5JY - ~TT for N f = 2 ■ ( 5 - 3 ) 



A^ 
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Specializing the general formula for the critical cou- 
pling a CTi R from the one-gluon exchange approximation 
to the Dyson- Schwinger equation, Eq. (|9.3j) (see Ap- 
pendix II) for the present case where R is the adjoint 
representation, one obtains a cr ^ a dj = n/(3N). Formally 
setting a>iR,2£,adj = a> a dj,cr yields the corresponding esti- 
mate for the critical number iV/ )Cr = 83/40 = 2.075. This 
may be rounded off to the nearest integer, giving Nf jCr for 
the adjoint representation. In view of the theoretical un- 
certainty in such an estimate, due to the strong-coupling 
nature of the physics involved, an SU(A r ) gauge theory 
with Nf = 2 adjoint fermions could be either slightly in- 
side the chirally broken, confined side of Nf >cr or slightly 
on the other side, where the theory is chirally symmet- 
ric and the evolution into the infrared is governed by an 
exact conformal IR fixed point. 

For the present case of Nf = 2 fermions in the adjoint 
representation of SU(JV), the coefficients of the beta func- 
tion are bi = N,b 2 = -107V 2 , b 3 = -(101/2)iV 3 , and 



6 4 = N 



1843 



-N - 312 



4((3)N 2 (N 2 + 72) . (5.4) 



At the four-loop level, the beta function has three ze- 
ros away from the origin, one of which is the four-loop IR 
zero, denoted otiR : 4£,adj- For N = 2, the others form an 
unphysical complex-conjugate pair, while for the other 
values of N that we consider, the others consist of a neg- 
ative one and a another, denoted uw,m which is not rel- 
evant to our study, since it is not reached by evolution of 
the coupling starting at small a for large [i. We list the 
numerical values of these zeros in Table IVHI 

The coefficients ci in Eq. (I2.5[) for 7 for this case are 
ca = 3A/(2tt), c 2 = (11A 2 )/(8tt 2 ), and c 3 = -N 3 /(2tt 3 ), 
with C4 given by 



4- 
7T C4 



N 



2 r 



192 



+ ^C(3)A 2 (iV 2 -9) . (5.5) 



(The term in C3 proportional to £(3) and the terms in 
C4 proportional to £(4) and £(5) vanish for the adjoint 
representation for arbitrary Nf.) 

Evaluating the two-loop expression in Eq. (|4.3I) for j m 
at the IR zero of the beta function, also calculated at the 
two- loop level, (Xm^adji we obtain 



l2t,adj{otIR,2l,adj) 



(11 - 4iV/)(277 - 34A/ + 40iV|) 



6(-17+ IQNff 



so that for the Nf — 2 case of interest here, 



(5.6) 



72e,adj{(XlR,2e,a.d 



41 

50 



= 0.820 for N f = 2 . (5.7) 



It is also of interest to evaluate the two-loop 7 m at 
the value of a cr from the one-gluon exchange (ladder) 
approximation to the Dyson-Schwinger equation. With 
N f = 2, this yields 



47 



(5.8) 



Evaluating the three-loop result for "f m at the IR zero 
of the beta function calculated at the three-loop level, 
ctiR,3i,adj, for the Nf = 2 case of interest, we obtain 



(5.9) 



which is again independent of N . At the four-loop level, 
the value of ju^djiaiiR^adj) does depend slightly on N. 
We list the values of these anomalous dimensions in Ta- 
ble IVIIII The most recent simulations of a lattice gauge 
theory with SU(2) gauge group and Nf — 2 fermions in 
the adjoint representation report 7 m = 0.49 ± 0.13 [23| . 
This is in agreement with the calculations of j m here at 
the three- and four-loop level, to within the uncertainties 
of the respective calculations. 



VI. SYMMETRIC AND ANTISYMMETRIC 
RANK- 2 TENSOR REPRESENTATIONS 

In this section we consider the SU(7V) theory with Nf 
fermions in the symmetric or antisymmetric rank-2 rep- 
resentation, denoted S2 and A2. Since a number of for- 
mulas are similar for these two cases, we will often give 
these in a unified way for both cases, denoted T2 (for 
rank-2 tensor representation), with ± signs distinguish- 
ing them. For S2, our analysis applies for any TV, while 
for A2, we restrict to N > 4, since the A2 representation 
is the singlet for SU(2) and is equivalent to the conju- 
gate fundamental representation for SU(3). Note that 
for SU(4), the A2 representation is self-conjugate. Also, 
since for SU(2) the S2 representation is the same as the 
adjoint representation, which has already been analyzed, 
we only consider the illustrative values N = 3, 4. 

For the two T2 cases, the general expression for the 
maximal value of Nf allowed by the requirement of 
asymptotic freedom, Eq. ()3.2|) . reduces to 



N 



UN 



f,max,T2 



2(N±2) ' 



(6.1) 



where the ± refers to S2 and A2, respectively. As N 
increases from 2 to oo, Nf tmaXj s2 increases monotonically 
from 2.75 to 11/2 = 4.5, and as N increases from 3 to oo, 
Nf,max.A2 decreases monotonically from 16.5 to the same 
limit, 4.5. The physical values of Nf^ max in both cases 
are the greatest integral parts of these rational numbers. 

For these representations, the general expression in Eq. 
(|3.5[) for the value of Nf at which the beta function co- 
efficient &2 changes sign from positive to negative with 
increasing Nf takes the form 



N 



17 TV 2 



f,b2z,T2 



(7V±2)(8iV±3-6iV- 1 ) 



(6.2) 



As a consequence of the general inequality (|3.6p , it follows 
that Nf yb2z ,s2 < Nf^ maXt s2 and Nf t i 2 z,A2 < Nf^ maXyA 2. 
For N — 2, the S2 representation is just the adjoint rep- 
resentation, so we only consider the illustrative values 
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N = 3, 4. The respective intervals Nf t b2z,S2 < Nf < 
Nf ima x,S2 for which the SU(A) gauge theory is asymptot- 
ically free and has an IR zero of f5 are 1.06 < Nf < 2.75 
for TV = 3 and 1.22 < Nf < 3.30 for N = 4. These 
ranges imply that the only physical integral values of Nf 
satisfying these conditions are Nf = 2, 3 for both SU(3) 
and SU(4). 

For large TV, Nf^2z,T2 has the series expansion 

17 323 6137 103547 n / 1 \ ln . 
Nf, b 2 Z<T2 = ¥ T ¥ ^ + ¥ ^T WW +0[ w ) (6.3) 

As TV increases from 2 to oo, Nf^2z,S2 increases mono- 
tonically from 17/16 = 1.0625 to 17/8 = 2.125, and as 
TV increases from 3 to oo, Nf y b2z,A2 decreases monoton- 
ically from 8.05 to the same limit, 2.125. This limit is 
twice the (TV-independent) value of Nf^zadj = 17/16 
for the adjoint representation. Thus, for large TV, the 
range (|3.7[) where the SU(TV) theory with Nf fermions in 
the S2 or A2 representation is asymptotically the same 
for both, namely, 17/8 < TV/ < 11/2; restricting TV/ to 



physical, integer values, this range consists of the three 
values TV/ = 3, 4, 5 

For our further discussion we assume that TV/ is in the 
range TV/ ib2 z,T2 < TV/ < Nf^ maXtT 2 where the theory is 
asymptotically free and the two-loop beta function has 
an IR zero, for the respective cases S2 and A2. This zero 
occurs at the value 

2tt(11A- 2A/(A±2) 
a m ,2e,T2 - _ 17N 2 + ^(gTV 2 ± 19TV =f 12TV" 1 ) ' 

(6.4) 

We have calculated am <n ^s2 and aiR iU i t A2 up to n = 
4 loops and list the results in Table IIXI and Table IXl 
The resultant /3 DS estimates for TV/ )Cr in the case of 
the S2 representation and N = 2, 3, 4 are A/ iCr .s2 = 
2.1, 2.5, 2.8, respectively. For the A2 representation 
with N = 4, one has Nf_ cr ^A2 — 8.1. 

The two-loop expression for the anomalous dimension, 
evaluated at a = aiR^i,T2, is 



l2l,T2{&IR,2l,T2) — 

(TV ± 2) (TV T 1) [UN - 2(TV ± 2)TV/ 



TV(554TV 2 ± 99TV - 198) + (-34TV 3 ± 22TV 2 =p 360) A/ + 20A(A ± 2) 2 A 2 
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- 17A 3 + (A ± 2)(8TV 2 ± 3A - 6) A/ 



(6.5) 



We list values of "f2t,S2{(XiR,2t,S2) for A = 2, 3, 4 in Table 
IXlland values of "f2i,A2{c t iR,2i,A2) for A = 4 in Table IXTfl 
with t = 2, 3. 

It is also of interest to evaluate the two-loop expression 
for 7 at the estimated a = a cr ^T2- This yields 



322A 2 ± 225TV - 450 - 10TV(TV ± 2)TV/ 

72,,T2K,,T2) 432(A±2)(A T 1) " 

(6.6) 

We list these values in Tables IXII and IXIII 

Evaluating the two-loop anomalous dimensions at the 
two-loop IR zero of the beta function, ^t^io-i R,2i,T2) , 
for TV/ equal to the respective (3 DS-estimated critical 
values, we obtain (again with T2 and the ± signs referring 
respectively to S2 and A2) 

72e,T2(otiR,2e,T2)\N f = N ft „. iT2 = 

374A 4 ± 578A 3 - 931TV 2 =f 900A + 900 



144(TV ± 2)(TV =f 1)(4A 2 ± 3A - 6) 
This has the large- TV expansion 

l2i,T2{ a IRae.,T2)\N f =N JiCT ,T2 = 



(6.7) 



187 

288 



17 n fJ_ 

128A + 1 A 2 



(6.8) 



The leading term has the value 187/288 ~ 0.649. 

From a lattice study of SU(3) gauge theory with Nf = 
2 fermions in the S2 (sextet) representation, Ref. [18| 
found that this theory is characterized by slow running 
behavior consistent with an (exact or approximate) IR 
fixed point, and further reported that j m < 0.6 where 
it was measured. For SU(3), the estimate of Q C t-,S2 in 
Eq. (j9.3p gives a cr ,s2 = 7r/10 = 0.31. Our results for 
the IR zero of /3 and the value of 7 m at this zero for 
TV = 3 and Nf = 2 are listed in Tables |IX] and El We 
find that aiR <n £ t s2 is approximately 0.84 at n — 2 loop 
level and decreases somewhat to 0.50 at three-loop level. 
The two-loop result for "f m is unphysically large, while 
the three-loop value of "f m at the corresponding three- 
loop IR zeros of j3 is about 1.3. These are somewhat 
larger than the values reported in Ref. [la ], although 
in assessing this comparison, one must take account of 
the significant strong-coupling uncertainties in our cal- 
culation stemming from the fact that otiR t s2 ~ 0(1). 
Our evaluation of the two-loop expression for j m at the 
ladder-Dyson-Schwinger estimate of a cr ^s2, is 0.65. 
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VII. EFFECTS OF NONZERO FERMION 
MASSES 

The global chiral symmetry that is operative if the 
fermions are massless, and the way that it is broken by 
fermion condensates, is well-known, and we do not review 
it here. However, it is worthwhile to comment on the sit- 
uation in which some fermion masses are nonzero. In this 
paper we generally assume that the fermions have zero in- 
trinsic masses in the Lagrangian describing the high-scale 
physics, and the only masses that they acquire arise dy- 
namically if they are involved in condensates that form as 
the gauge interaction becomes sufficiently strongly cou- 
pled in the infrared. This is a well-motivated assumption 
if the vectorial gauge theory arises as a low-energy effec- 
tive field theory from an ultraviolet completion which is 
a chiral gauge theory. In turn, this is natural if the latter 
theory becomes strongly coupled, since it can then form 
fermion condensates that self-break it down to the vecto- 
rial subgroup symmetry. However, one may also choose 
to focus on the vectorial gauge theory as an ultraviolet- 
complete theory in itself. In a vectorial gauge theory, an 
intrinsic (bare) mass term for a fermion ip, C m = —rmj)'^, 
is allowed by the gauge invariance. Hence, one may con- 
sider a more general situation in which the fermions may 
have such intrinsic (hard) masses in the high-scale La- 
grangian [44] . In this the reference scale /x de- 
creases below the value of the hard mass of some fermion 
m,f, the beta function changes from one that includes this 
to one that excludes this fermion. If the hard fermion 
masses are small compared with the scale A in the situ- 
ation where the theory confines and breaks chiral sym- 
metry spontaneously, then these hard masses have only a 
small effect. However, if some of the hard fermion masses 
are sufficiently large, then as \i decreases below their scale 
and the corresponding fermions are integrated out of the 
low-energy theory below this scale, this can significantly 
change the infrared properties of the resultant theory. 

In applications of slowly running gauge theories 
to technicolor theories, at the scale Axe where the 
SXJ(Ntc) gauge coupling grows to 0(1) and is influenced 
by the presence of an approximate IR zero of the TC 
beta function, there can also be non-negligible effects 
due to four-fermion operators arising from the higher- 
lying extended technicolor dynamics [J]-@; fToj|- |45l. 14^. 
and these can affect the scaling properties of -tptp. Similar 
comments apply for topcolor-assisted technicolor 0, S3 • 

VIII. CONCLUSIONS 

In this paper we have studied the evolution of an 
asymptotically free vectorial SU(JV) gauge theory from 



high scales to the infrared taking account of higher-loop 
corrections to the beta function and the anomalous di- 
mension 7 m for fermions in the fundamental, adjoint, and 
rank-2 symmetric and antisymmetric representations S2 
and A2. We have compared our results with lower-order 
calculations. We have shown that, for fixed N and Nf, 
in the range for which the two-loop beta function has an 
IR zero, the value of this zero decreases as one goes from 
the two-loop to the three-loop calculations, and we have 
determined this decrease quantitatively. Going further, 
we have shown that there is a smaller fractional increase 
in the value of this IR zero when calculated to four-loop 
accuracy, with the final four-loop result still smaller than 
the two-loop value. We have analyzed instanton effects 
and have demonstrated that they tend to increase the 
value of the IR zero of the beta function somewhat. A 
major part of our work has been the evaluation of the 
anomalous dimension 7 m of iptp at the IR zero of the beta 
function at the £ = 2, 3, 4 loop levels. This zero is approx- 
imate or exact, depending on whether for a given N, the 
value of Nf is below or above the critical value iV/ jCr be- 
low which there is spontaneous chiral symmetry breaking 
associated with the formation of a fermion condensate. 
We have found that this 7 m at the (approximate or ex- 
act) IR zero of the beta function decreases as one goes 
from two-loop to three-loop order, and that the four-loop 
values also tend to be somewhat less than those at the 
two-loop level. The values that we have calculated for 7 m 
at the IR zero of the beta function tend to be somewhat 
smaller than unity. We have compared our higher-loop 
calculations with results from recent lattice simulations 
and have found general agreement. We believe that the 
higher-loop calculations reported here should provide a 
useful reference for comparison with ongoing and future 
lattice measurements. 
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IX. APPENDIX I 



For the reader's convenience, we list the three-loop 
beta function coefficient, in the MS scheme [Hj], 



13 



h = + +T jNf 



2Cj 



205 



C A C f 



1415 
27 



(TfN f ) 2 



44 



-C t 



158 
~27 



C A 



(9.1) 



The four-loop coefficient is given in Ref. 29] and is a cubic polynomial in Nj. We note that the coefficients of the 
N® (which is independent of the fermion representation) is positive, and the coefficient of the Nj term is positive for 
an arbitrary fermion representation. 

Our normalizations for the quadratic Casimir and trace invariants of a Lie group are standard. The quadratic 
Casimir invariant Ca(-R) for the representation R is given by X)a=i X^Hi [Dn(T a )]ij[Dii(T a )]jk — C2{R)8ik, 
where a, b are group indices, o(G) is the order of the group, T a are the generators of the associated Lie alge- 
bra, and Dn(T a ) is the matrix form of the T a in the representation R. The trace invariant T(R) is defined by 

T^WRiTaMDniT^ii = T(R)8 ab . 

From the calculations of the coefficients of the perturbative expansion of the anomalous dimension 7 m in the MS 
scheme to four- loop order in Ref. [3(j, we record the three-loop coefficient 



C3 



2C f 



129 r2 



129 

—CfC A 



11413 
108 



Ci + C / T / 7V / (-46 + 48C(3)) 



C A T f N f ( 5 -§ 



(9.2) 



We have used the four-loop coefficient C4 from Ref. [3C 



for our calculations, but it is too lengthy to reproduce here. 
I 



A. Appendix II: Beta-Dyson-Schwinger Estimate 

Of Nf >or 

In this appendix we briefly review the ft DS estimate 
of Nf >cr . In the one-gluon exchange (also called ladder) 
approximation to the Dyson- Schwinger equation for the 
fermion propagator with an initially massless fermion in 
the representation R of the gauge group, one finds a solu- 
tion with a dynamically generated, nonzero fermion mass 
if the coupling ol{h) exceeds a critical value a cr r given 
byGLIEl 



a cr ,R 



7T 

3C 



/ 



(9-3) 



Setting this equal to the two-loop expression for the IR 
zero of ft then yields an estimate for Nf lCr to this order, 
namely 



Nt = 

iv f,cr — 



C A (66C f + 17C A ) 
10T f (C A +3C f ) 



(9.4) 



We call this the ftDS estimate of Nf cr since it combines a 
calculation of am from the ft function with the estimate 
of a crt R from the ladder approximation to the Dyson- 
Schwinger equation for the fermion propagator. In the 
same ladder approximation, one finds 7 m = 1 at a = 
ot cr ,R [J (which also holds for the DS analysis at a UV- 
stable fixed point (48[). For the gauge group SU(iV) with 
the illustrative values of N used for the tables, namely 
N = 2, 3, 4, Nf^ cr j un d is equal to 7.9, 11.9, and 15.9, 
respectively, with the large-TV form Nf ~ AN. For S2, 
the symmetric rank-2 tensor representation, N = 2, 3, 4, 
Nf yCr ,s2 is equal to 2.075, 2.5, and 2.9, increasing toward 
the limit 11/2 = 5.5 in the large- N limit. In the case of 
A2, the antisymmetric rank-2 tensor reprepresentation, 



for N = 3, the result is the same as for the fundamental 
representation, while for N = 4, one has Nf >cr it t A2 — 8.1, 
and as N — > 00, Nf lCr u,A2 decreases toward the limit 
11/2. 

One understands that, a priori, there could be sig- 
nificant uncertainty in these estimates because of the 
strong-coupling nature of the physics involved and the 
one-gluon approximation used for the solution of the 
Dyson- Schwinger equation. Moreover, the DS equation 
analysis is semi-perturbative in the sense that it contains 
polynomial dependence on a, and it neglects nonpertur- 
bative effects associated with confinement and instan- 
tons. However, corrections to the one-gluon exchange 
approximation have been analyzed and found not to be 
too large 0. Recent lattice simulations for SU(3) are 
in broad agreement, to within the uncertainties, with the 
above prediction of N fyCr ~ 12 HH-dl, 0]. Some of the 
success of the ftDS prediction for Nf_ cr may arise from 
the fact that two major physical effects that it ignores, 
namely confinement and instantons, would shift Nf jCr in 
opposite directions and hence tend to cancel each other 
out M. 



X. APPENDIX III - PADE RESULTS 

In this appendix we collect some relevant results on 
Pade approximants. Given a Taylor (or asymptotic) se- 
ries expansion around z = for the function f(z), 



f(z) = £ /„z" 



(10.1) 



n=0 



one can construct a set of [p, q] Pade approximants, 
namely rational functions comprised of a numerator poly- 
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nomial of degree p and a denominator polynomial of 
degree q, such that p + q = n max — 1, of the form 
=oPi zJ )/(Sfc=o9fc zfc )- Without loss of generality, 
one can divide numerator and denominator by qo, so that, 
after redefinition of the coefficients, one has 

m ' <z) = ttES^- <ia2) 



The p + q + 1 coefficients pj with < j < p and q^ with 
1 < /c < q are uniquely determined in terms of the /„ 
coefficients with < n < n max by expanding the [p, q] 
Pade approximant in a Taylor series around z = and 
solving the set of n max linear equations. 
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TABLE I: Values of the £-loop beta function coefficients bg defined 
in Eq. J2.2B in the SXJ(N) gauge theory with Ne fermions trans- 
forming according to the fundamental representation, as functions 
of N and Nf , for the range 1 13, II) where the theory is asymptotically 
free. 



N 


AT/ 


h 




bz 


h 


2 





0.584 


0.287 


0.213 


0.268 


2 


1 


0.5305 


0.235 


0.154 


0.191 


2 


2 


0.477 


0.184 


0.099 


0.127 


2 


3 


0.424 


0.132 


0.047 


0.078 


2 


4 


0.371 


0.080 


-0.0003 


0.044 


2 


5 


0.318 


0.0285 


-0.044 


0.024 


2 


6 


0.265 


-0.023 


-0.084 


0.020 


2 


7 


0.212 


-0.075 


-0.120 


0.030 


2 


8 


0.159 


-0.127 


-0.152 


0.057 


2 


9 


0.106 


-0.178 


-0.180 


0.099 


2 


10 


0.053 


-0.230 


-0.205 


0.156 


3 





0.875 


0.646 


0.720 


1.173 


3 


1 


0.822 


0.566 


0.582 


0.910 


3 


2 


0.769 


0.485 


0.450 


0.681 


3 


3 


0.716 


0.405 


0.324 


0.485 


3 


4 


0.663 


0.325 


0.205 


0.322 


3 


5 


0.610 


0.245 


0.091 


0.194 


3 


6 


0.557 


0.165 


-0.016 


0.099 


3 


7 


0.504 


0.084 


-0.118 


0.039 


3 


8 


0.451 


0.004 


-0.213 


0.015 


3 


9 


0.398 


-0.076 


-0.303 


0.025 


3 


10 


0.345 


-0.156 


-0.386 


0.072 


3 


11 


0.292 


-0.236 


-0.463 


0.154 


3 


12 


0.239 


-0.317 


-0.534 


0.273 


3 


13 


0.186 


-0.397 


-0.599 


0.429 


3 


14 


0.133 


-0.477 


-0.658 


0.622 


3 


15 


0.080 


-0.557 


-0.711 


0.852 


3 


16 


0.0265 


-0.637 


-0.758 


1.121 


4 





1.17 


1.15 


1.71 


3.50 


4 


1 


1.11 


1.04 


1.46 


2.88 


4 


2 


1.06 


0.932 


1.22 


2.31 


4 


3 


1.01 


0.824 


0.986 


1.80 


4 


4 


0.955 


0.716 


0.762 


1.36 


4 


5 


0.902 


0.607 


0.546 


0.972 


4 


6 


0.849 


0.499 


0.339 


0.647 


4 


7 


0.796 


0.391 


0.140 


0.385 


4 


8 


0.743 


0.283 


-0.051 


0.184 


4 


9 


0.690 


0.175 


-0.234 


0.046 


4 


10 


0.637 


0.066 


-0.409 


-0.029 


4 


11 


0.584 


-0.042 


-0.575 


-0.040 


4 


12 


0.531 


-0.150 


-0.733 


0.013 


4 


13 


0.477 


-0.258 


-0.883 


0.131 


4 


14 


0.424 


-0.366 


-1.025 


0.314 


4 


15 


0.371 


-0.474 


-1.16 


0.562 


4 


16 


0.318 


-0.583 


-1.28 


0.877 


4 


17 


0.265 


-0.691 


-1.40 


1.26 


4 


18 


0.212 


-0.799 


-1.51 


1.71 


4 


19 


0.159 


-0.907 


-1.61 


2.22 


4 


20 


0.106 


-1.015 


-1.70 


2.81 


4 


21 


0.053 


-1.124 


-1.79 


3.46 
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TABLE II: Values of N f>b2z , Nf, b3z> ±, and N fMzJ , i = 2,3, 
for SU(iV) with Nf fermions in the fundamental representation. 
We only list physical, i.e., real, non-negative values. Thus, since 
Nf,bz4.i < 0, is is not included. 



N 


^ f,max 


Nf,b2z 


(Nj\fiz,- , Nf,b3z,+ ) 


(Nf t b4z,2, Nf,b4z,3) 


2 


ll 


5.55 


(3.99, 27.6) 


none 


3 


16.5 


8.05 


(5.84, 40.6) 


none 


4 


22 


10.61 


(7.73, 53.8) 


(9.51,11.83) 



TABLE III: Values of the (approximate or exact) IR zeros in a 
of the SU(iV) beta function with Nf fermions in the fundamental 
representation, for N = 2, 3, 4, calculated at ra-loop order, and 
denoted as cijr n £. For each N, we only give results for the integral 
Nf values in the range 113. 7t where the theory is asymptotically 
free and the two-loop beta function has an infrared zero. For the 
four-loop beta function, the cubic equation H3.32H has three zeros, 
one of which is negative, one of which is ckjr 41, and the third of 
which is positive but farther from the origin. We include the latter, 
denoted as u . We also list zeros from the [1,2] and [2,1] Pade 
approximants to the four-loop beta function. 



V 


N f 


OlIR,2t 


CtlR,3l 




Q /iJ,«,[l,2] 


&IRAl,[2,\\ 


Ct4l, u 


2 


7 


2.83 


1.05 


1.21 


2.30 


1.16 


4.12 


2 


8 


1.26 


0.688 


0.760 


0.952 


0.741 


3.11 


2 


9 


0.595 


0.418 


0.444 


0.475 


0.438 


2.395 


2 


10 


0.231 


0.196 


0.200 


0.202 


0.200 


1.97 


3 


10 


2.21 


0.764 


0.815 


1.47 


0.807 


5.62 


3 


11 


1.23 


0.578 


0.626 


0.871 


0.616 


3.29 


3 


12 


0.754 


0.435 


0.470 


0.561 


0.462 


2.295 


3 


13 


0.468 


0.317 


0.337 


0.367 


0.333 


1.78 


3 


14 


0.278 


0.215 


0.224 


0.231 


0.222 


1.48 


3 


15 


0.143 


0.123 


0.126 


0.127 


0.125 


1.29 


3 


16 


0.0416 


0.0397 


0.0398 


0.0398 


0.0398 


1.15 


4 


13 


1.85 


0.604 


0.628 


1.14 


0.625 


6.94 


4 


14 


1.16 


0.489 


0.521 


0.776 


0.516 


3.49 


4 


15 


0.783 


0.397 


0.428 


0.556 


0.422 


2.30 


4 


16 


0.546 


0.320 


0.345 


0.407 


0.340 


1.73 


4 


17 


0.384 


0.254 


0.271 


0.298 


0.267 


1.40 


4 


18 


0.266 


0.194 


0.205 


0.215 


0.203 


1.19 


4 


19 


0.175 


0.140 


0.145 


0.149 


0.145 


1.05 


4 


20 


0.105 


0.091 


0.092 


0.0930 


0.0921 


0.947 


4 


21 


0.0472 


0.044 


0.044 


0.0444 


0.0443 


0.870 
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TABLE IV: Estimates of a C r,R from the one-gluon exchange ap- 
proximation to the Dyson-Schwingcr equation for the fermion prop- 
agator. Values are listed for SU(JV) with 2 < N < 6 and the rep- 
resentations R = (i) fundamental (fund), (ii) adjoint (adj), (iii) 
symmetric rank-2 tensor (S2), and (iv) antisymmetric rank-2 (A2). 



N 


OCcr,fund 


a cr ,adj 


OLcr,S2 


Otcr,A2 


2 


1.40 


0.52 


0.52 




3 


0.79 


0.35 


0.31 


0.79 


4 


0.56 


0.26 


0.23 


0.42 


5 


0.44 


0.21 


0.19 


0.29 


6 


0.36 


0.17 


0.16 


0.22 
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TABLE V: Values of the £-loop coefficients in the series expan- 
sion 1 12. 51 1 for the anomalous dimension j m , as functions of N and 
Nf, for the range 1 13.111 where the theory is asymptotically free. 



N 




Cl 


C2 


C3 


C4 


2 





0.358 


0.318 


0.310 


0.329 


2 


1 


0.358 


0.302 


0.254 


0.234 


2 


2 


0.358 


0.286 


0.195 


0.143 


2 


3 


0.358 


0.270 


0.134 


0.0577 


2 


4 


0.358 


0.254 


0.0712 


-0.0218 


2 


5 


0.358 


0.239 


0.00656 


-0.0952 


2 


6 


0.358 


0.223 


-0.0601 


-0.162 


2 


7 


0.358 


0.207 


-0.129 


-0.222 


2 


8 


0.358 


0.191 


-0.199 


-0.274 


2 


9 


0.358 


0.175 


-0.272 


-0.319 


2 


10 


0.358 


0.1595 


-0.346 


-0.355 


3 





0.637 


0.853 


1.26 


2.03 


3 


1 


0.637 


0.825 


1.11 


1.64 


3 


2 


0.637 


0.796 


0.957 


1.27 


3 


3 


0.637 


0.768 


0.801 


0.909 


3 


4 


0.637 


0.740 


0.642 


0.561 


3 


5 


0.637 


0.712 


0.479 


0.227 


3 


6 


0.637 


0.684 


0.312 


-0.0926 


3 


7 


0.637 


0.656 


0.142 


-0.396 


3 


8 


0.637 


0.628 


-0.0313 


-0.683 


3 


9 


0.637 


0.599 


-0.208 


-0.953 


3 


10 


0.637 


0.571 


-0.389 


-1.21 


3 


11 


0.637 


0.543 


-0.573 


-1.44 


3 


12 


0.637 


0.515 


-0.760 


-1.65 


3 


13 


0.637 


0.487 


-0.951 


-1.85 


3 


14 


0.637 


0.459 


-1.145 


-2.02 


3 


15 


0.637 


0.431 


-1.34 


-2.18 


3 


16 


0.637 


0.402 


-1.54 


-2.31 


4 





0.895 


1.60 


3.17 


6.86 


4 


1 


0.895 


1.56 


2.89 


5.88 


4 


2 


0.895 


1.52 


2.61 


4.93 


4 


3 


0.895 


1.48 


2.32 


4.00 


4 


4 


0.895 


1.44 


2.03 


3.09 


4 


5 


0.895 


1.40 


1.73 


2.21 


4 


6 


0.895 


1.36 


1.43 


1.36 


4 


7 


0.895 


1.33 


1.12 


0.526 


4 


8 


0.895 


1.29 


0.808 


-0.275 


4 


9 


0.895 


1.25 


0.492 


-1.05 


4 


10 


0.895 


1.21 


0.170 


-1.79 


4 


11 


0.895 


1.17 


-0.157 


-2.50 


4 


12 


0.895 


1.13 


-0.488 


-3.18 


4 


13 


0.895 


1.09 


-0.825 


-3.83 


4 


14 


0.895 


1.05 


-1.17 


-4.45 


4 


15 


0.895 


1.01 


-1.51 


-5.025 


4 


16 


0.895 


0.969 


-1.86 


-5.57 


4 


17 


0.895 


0.930 


-2.22 


-6.08 


4 


18 


0.895 


0.890 


-2.58 


-6.54 


4 


19 


0.895 


0.850 


-2.95 


-6.97 


4 


20 


0.895 


0.811 


-3.32 


-7.37 


4 


21 


0.895 


0.771 


-3.69 


-7.71 



20 



TABLE VI: Values of the anomalous dimension in the SU(AT) 
theory with Nf fermions in the fundamental representation, y m , 
calculated to the n-loop order in perturbation theory and evaluated 
at the IR zero of the beta function calculated to this order, cem n £, 
for i = 2,3,4. We denote these as ■y n e(aiR,ni)- For sufficiently 
small Nf > Nf b 2z m each N case, Oijji ^l ls so large that the 
formal value of f2l( Q lR it) is larger than 2 and hence unphysical; 
we indicate this by placing these values in parentheses. 



N 


JV / 




^f3t\ a IR,3£ ) 


y4£\QlR,4£} 


2 


7 


(2 67) 


0.457 


0325 


2 


8 


0.752 


0.272 


0.204 


2 


9 


0.275 


0.161 


0.157 


2 


10 


0.0910 


0.0738 


0.0748 


q 
o 


i n 

-LU 


f4 1 Q"! 

V / 


D R47 


U. LOU 


q 


]^ 


1 fil 


fl 4^Q 


n 9^0 


q 


1 9 


u. it o 


n 9 


U.i JO 


q 
o 




n 4fl4 


n 99n 


n 91 n 


3 


14 


0.212 


0.146 


0.147 


3 


15 


0.0997 


0.0826 


0.0836 


3 


16 


0.0272 


0.0258 


0.0259 


4 


13 


(5.38) 


0.755 


0.192 


4 


14 


(2.45) 


0.552 


0.259 


4 


15 


1.32 


0.420 


0.281 


4 


16 


0.778 


0.325 


0.269 


4 


17 


0.481 


0.251 


0.234 


4 


18 


0.301 


0.189 


0.187 


4 


19 


0.183 


0.134 


0.136 


4 


20 


0.102 


0.0854 


0.0865 


4 


21 


0.0440 


0.0407 


0.0409 



TABLE VII: Values of the (approximate or exact) IR zeros in a 
of the SXJ(N) beta function with Nt = 2 fermions in the adjoint 
representation, for N = 2,3, 4, calculated at n-loop order, and de- 
noted as ciiR,nl,adi- F° r the four-loop beta function, the cubic 
equation l|3.32[l has three zeros, one of which is is aiR,4t, a dj ■ De- 
pending on N, there may be another real zero, denoted om u a ^j, 
at a larger value of a. We also list zeros from the [1,2] and [2,1] 
Pade approximants to the four-loop beta function. 



N 


OLIR,2t,adj 


OLIR,Zt,adj 


dIR,4l,adi 


aiR,4e,[l,2],adj 


&IR,U,\2,l},a,dj 


0!4e,u,adj 


2 


0.628 


0.459 


0.450 


0.455 


0.449 




3 


0.419 


0.306 


0.308 


0.317 


0.308 


9.38 


4 


0.314 


0.2295 


0.234 


0.242 


0.233 


3.29 



21 



TABLE VIII: Values of the anomalous dimension 7 m in an SU(iV) 
gauge theory with Nf = 2 (Dirac) fermions in the adjoint repre- 
sentation, calculated to the n-loop order in perturbation theory 
and evaluated at the IR zero of the beta function calculated to 
this order, for n = 2,3,4. We denote these as "Yne^dji&l R,nl,adj)- 
We also list the value of ^e, adj evaluated at a equal to the /3DS 
estimate, Eq. 1 19.311 , for a criac ij). 



N 


l2l,adj{0!lRfil,adj) 


y3e,adj(0!lR^e,adj) 


y4e,adj{ctIR,4t,adj) 


72f,adj(Q! cr , a[ Jj) 


2 


0.820 


0.543 


0.500 


0.653 


3 


0.820 


0.543 


0.523 


0.653 


4 


0.820 


0.543 


0.532 


0.653 



TABLE IX: Values of the (approximate or exact) IR zero in a of 
the SXJ(N) beta function with Nf = 2 fermions in the symmetric 
rank-2 (i.e., S2) representation, for N = 3,4, calculated at n-loop 
order, and denoted as ct IR n i g2- 



N 


N f 


OtIR,2l,S2 


aiR,3t,S2 


CtIRM,S2 


3 


2 


0.842 


0.500 


0.470 


3 


3 


0.085 


0.079 


0.079 


4 


2 


0.967 


0.485 


0.440 


4 


3 


0.152 


0.129 


0.131 



TABLE X: Values of the (approximate or exact) IR zero in a of 
the SU(4) beta function with Nf fermions in the antisymmetric 
rank-2 (i.e., A2) representation, for the range 5 < Nf < 10 where 
the theory is asymptotically free and has an IR zero of the beta 
function, calculated at n-loop order, and denoted as otiR in t t A2- 



N 


N f 


OiIR,2l,A2 


OiIR,Zl,A2 


OiIRAl,A2 


4 


6 


2.17 


0.664 


0.770 


4 


7 


0.890 


0.437 


0.502 


4 


8 


0.449 


0.287 


0.319 


4 


9 


0.225 


0.174 


0.184 


4 


10 


0.090 


0.080 


0.082 



TABLE XI: Values of 7 m in an SV(N) gauge theory with Nf 
fermions in the symmetric rank-2 tensor representation S2, cal- 
culated to the n-loop order in perturbation theory and evaluated 
at the IR zero of the beta function calculated to this order, for 
n = 2,3,4. We denote these as ■y„i S2( a lR nl,S2)- We also list 
72£, S2 evaluated at a equal to the estimate Eq. I|9.3[l for ct cr ,S2- 



N 


N f 


l2t,S2{(XlR,2t,S2) 


l3l,S2{&IR,3,l,S2) 


7«,S2(«/fl,4^,S2) 


72£,S2(«cr,S2) 


3 


2 


(2.44) 


1.28 


1.12 


0.653 


3 


3 


0.144 


0.133 


0.133 


0.619 


4 


2 


(4.82) 


(2.08) 


1.79 


0.659 


4 


3 


0.381 


0.313 


0.315 


0.629 



22 



TABLE XII: Values of r y m in an SU(JV) gauge theory with Nf 
fermions in the antisymmetric rank-2 tensor representation A2, cal- 
culated to the n-loop order in perturbation theory and evaluated 
at the IR zero of the beta function calculated to this order, for 
N = 4 and n = 2,3,4. We denote these as ■y n i,A2(c*l R,nl,A2)- We 
also list 72£,A2 evaluated at a equal to the estimate Eq. 119, 3H for 
a cr ,A2- 



N 


Nf 


J2£.,A2(CXIR,2£,A2) 


~l3l,A2{a-IR,3l,A2) 


r Y4,l,A2\OllR,4l,A2) 


72<?,A2(a C r,A2) 


4 


6 


(9.78) 


1.38 


0.293 


0.769 


4 


7 


(2.19) 


0.695 


0.435 


0.750 


4 


8 


0.802 


0.402 


0.368 


0.732 


4 


9 


0.331 


0.228 


0.232 


0.713 


4 


10 


0.117 


0.101 


0.103 


0.695 



